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HEAT  TRANSFER  IN  ROCKET  ENGINE 
COMBUSTION  CHAMBERS  AND  NOZZLES 


P.  G.  Anderson,  G.  C.  Cheng  and  R.  C.  Fanner 
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Abstract 


Complexities  of  liquid  rocket  engine  heat  transfer  which  involve  the  injector  faceplate  and 
regeneratively  and  film  cooled  walls  are  being  investigated  by  computational  analysis.  A 
conjugate  heat  transfer  analysis  will  be  used  to  describe  localized  heating  phenomena  associated 
with  particular  injector  configurations  and  coolant  channels  and  film  coolant  dumps.  These 
components  are  being  analyzed,  and  the  analyses  verified  with  appropriate  test  data.  Finally 
the  component  analyses  will  be  synthesized  into  an  overall  flowfield/heat  transfer  model.  The 
FDNS  code  is  being  used  to  make  the  component  analyses.  Particular  attention  is  being  given 
to  the  representation  of  the  thermodynamic  properties  of  the  fluid  streams  and  to  the  method  of 
combining  the  detailed  models  to  represent  overall  heating.  Unit  flow  models  of  specific  coaxial 
injector  elements  have  been  developed  and  will  be  described. 


Since  tdst  data  from  the  NLS  development  program  are  not  available,  new  validation  heat 
transfer  data  has  been  sought.  Suitable  data  was  obtained  from  a Rocketdyne  test  program  on 
a model  hydrocarbon/oxygen  engine.  Simulations  of  this  test  data  will  be  presented. 

Recent  interest  in  the  hybrid  motor  have  established  the  need  for  analyses  of  ablating 
solid  fuels  in  the  combustion  chamber.  Analysis  of  a simplified  hybrid  motor  will  also  be 
presented. 
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AXISYMMETRIC 
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/ pv»  At,  H„  are  functions  of  t and  p is  a function  of  t and 


MAIN  INJECTOR  ELEMENT  (PARTIAL  GRID) 
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MODIFICATIONS  TO  FDNS 
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(ITERATIVE  PROCEDURE  HAS 
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(3)  Uniform  O/F  Distribution 

( m,  = 4.45  Ib/sec,  rh2  = 13.35  Ib/sec,  m3  = 26.7  Ib/sec  ) 
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Grid  for  the  Axisymmetric  Configuration 
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Axisymmetric  Hybrid  Configuration 


Temperature  ,n  the  fuel  port  region  (degree  K) 
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PARAMETERS  AFFECTING  RADIATION 
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Tntprnal  Surface  Characteristics 

• Not  well  known 

• Assume  reflectivity  of  0.2 
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GAS  EMISSIVITY. 
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Wide  band  model  correlation  parameters  for  various  gases 
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PHASE  FUNCTION  AS 
DETERMINED  BY 
MIE  THEORY 
CALCULATION 


A120,  RADIATION  PROPERTIES 
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Original  A1205  Optical  Properties  in  SIRRM 
Original  A12Oj  Optical  Property  Data  With  MIE  Code 
OPTROCK  Properties  in  SIRRM 

OPTROCK  Properties  Calculated  With  A MIE  Code 
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Plume  Radiation  from  a Small  Scale  F-l  Engine 

(Reardon,  Radiation  Handbook) 
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FOR  SRM 

(from  Netzer’s  experiments) 
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PRESENT  PLUME  RADIATION  METHODOLOGY 
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the  YIX  method  (Tan  & Howell) 

• • uses  piecewise-constant  interpolation  and 
stores  intermediate  integral  evaluations 
to  improve  the  PIM  method. 
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Limitation:  For  N=1  or  3,  accurate  for  only  optically  thick  media. 


THE  MODIFIED  DIFFERENTIAL  APPROXIMATION 

The  MDA  generalizes  P,  solution  to  treat  arbitrary  optical  thickness. 
Let  surface  & media  contribute  to  intensity: 
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THE  IMPROVED  DIFFERENTIAL  APPROXIMATION 


1017 


o 

z 

HH 

C/3 

P 


c 

*3 

*-> 

O 

o 


00 

r- 

00 

s 


0) 

1 

P 

G 

cd 


00 


O 

J> 

O 

00 


ai 

u- 

<D 

Cu 


00 


C 

P 

cd 


<u 

o 

tcd 

u. 

P 

C/3 


<D 

w 

o 


E 

3 

00 

4-> 

00 

p 

£ 


Jp 

O 

15 

£ 


GO 

* 


<D 

£ 

cd 

SP 

<u 

c 

JD 

>> 

Jp 

w 

o 


8 j 

” Pi 

Q £ 

<N  £ 

cd  jd 

k O- 
.O  cd 

T3  w 
<L>  *_, 

t:  c 

°o 

S'  °* 

U -Q 

c *C 

<U  CxQ 

<u  17 

•°  -s 

SS  « 

2 o 

1 2 

.1  § 

*- > . — i 

— 3 

O u 

c/3  ^ 

•o  2 

^ -a 

cx  o 

3 So 
O 3 

« <D 

§ S3 

•S  * 

i£ 

o o 

h4 

S2 

.2  "5 

•S2 

ctf  £ 


1018 


CONCLUSIONS 


G 

<D 

<D 

& 

O IS 

^ to 
u on 
O 

Oh  O 


03  11 

' O yj 


m 

O £ 

3. 

<c  00 

Ui  vo 

O <2o 

8 

tl| 

8^ 
g fe 

<2  tn  O 

r-  cd 

^ O II 

(D  ^ 11 

"O  o vu 

O • 1-H 

C Ox  <D 
G o 
£ G 

TD  -p 
i a o jp 
C/D  D 

5 x>  •£  •£ 

O ;>  CVj  £ 

o I |g  £ 

6 £ S ^ 

— S3  g c 

5 G Gh  -S3 
H X3 

1)  <D  p 

G 3 £ 

* 


M 7? 

t_  <0 

§ & 
§ <D 

* 3 

C3 

5 § 

y 3 

S3  p 


1019 
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TZU-YI  YU , Bharat  K.  Soni  3 

NSF/ERC  For  Computational  Field  Simulation  , qj  j j ^ (~j  Cj 

Mississippi  State  University.  MS39762  ' 

ABSTRACT 

L.  fie  field  of  Grid  Generation  and  fie  CAD/CAM  . to  are  numerous  geomedy  fonMB  wUch 
the  designer  to  spend  a (pea,  deal  of  to  manipulating  geometrical  entities  in  onler  to  achieve  a 
nsefhi  sculptured  geometrical  descripdou  for  grid  generation  A,  so  in  fils  process,  to  is  a danger  of  loosing 
fideh^,  of  fie  geomedy  under  consideradon.  This  stresses  fie  impedance  of.  standard  geomedy  definition  for 
fie  communication  link  between  varing  CAD/CAM  and  grid  system  . The  ICES  (Initial  Graphics  Exchttnge 
Specification)  (Refl)  file  is  a widely  used  communication  between  CAD/CAM  and  the  analysis  tools  The  scien- 
tists ..  NASA  Resetuch  Centem-  including  NASA  AmesNASA  Inmgley.NASA  Lewis  and  NASA  Marshall 
- have  recognised  this  importance  and  therefore,  in  1992  they  fomted  the  committee  of  fie  'NASA-IGES"  which 
fie  subset  of  the  standard  IGES  This  committee  stresses  the  importance  and  encourage  the  CFD  community 
to  use  the  standani  IGES  file  for  fie  interface  between  the  CAD/CAM  and  CFD  analysis  . Also  two  of  the 
IGES  entities  - fie  NURBS  Crnve  (Entity  126)  and  NURBS  Surface  (Entity  128)  - which  have  many  useful 
geometric  properties  — like  the  convex  hull  property  . local  control  property  and  affine  invariance,  also  widely 
utilized  analytical  geometries  can  be  accurately  represented  using  NURBS.  This  is  improtant  in  today  grid  genera- 
tion  tools  because  of  the  emphasis  of  the  interactive  design. 

To  satisfy  the  geomedy  dansfonnation between  fie  CAD/CAMsystemandGridGeneralionfield  theCAGI- 
Compnter  Aided  Geomedy  Design  is  developed  . which  include  fie  Geomedy  Transfonnation  . Geomedy 

Manipulation  rod  Geomedy  Generates  well  as  fie  user  interface.  A self  explanatoypictoritivIewsofCAG. 
modules  and  links  is  shown  in  Figure  1 . 

This  paper  will  present  the  successful  development  IGES  file  tiansfonner  rod  application  of  NURBS 
definition  (Ref  3)  in  the  grid  generation  (Ref  4,5) 
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Geometry  Generation : 
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Abstract 

Many  CFD  grid  generation  systems  are  in  use  nationally,  but  few  comparative  studies 
have  been  performed  to  quantify  their  relative  merits.  A study  was  undertaken  to 
systematically  evaluate  and  select  the  best  CFD  grid  generation  codes  available. 
Detailed  evaluation  criteria  were  established  as  the  basis  for  the  evaluation 
conducted.  Descriptions  of  thirty  four  separate  criteria,  grouped  into  eight  general 
categones  are  provided.  Benchmark  test  cases,  developed  to  test  basic  features  of 
selected  codes,  are  described  in  detail.  Scoring  guidelines  were  generated  to 
establish  standards  for  measuring  code  capabilities,  ensure  uniformity  of  ratings,  and 
minimize  personal  bias  among  the  three  code  evaluators.  Ten  candidate  codes  were 
identified  from  government,  industry,  universities,  and  commercial  software 
companies.  A three  phase  evaluation  was  conducted.  In  Phase  1,  ten  codes 
identified  were  screened  through  conversations  with  code  authors  and  other  industry 
experts.  Seven  codes  were  carried  forward  into  a Phase  2 evaluation  in  which  all 
codes  were  scored  according  to  the  predefined  criteria.  Two  codes  emerged  as  being 
significantly  better  than  the  others;  RAGGS  and  GRIDGEN.  Finally,  these  two  codes 
were  carried  forward  into  a Phase  3 evaluation  in  which  complex  3-D  multizone  grids 
were  generated  to  verify  capability. 
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PHASE  2 - PRELIMINARY  EVALUATION 
PHASE  3 - FINAL  EVALUATION 
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INTERNAL  GEOMETRY  CREATION  CAPABILITY  (5/10) 
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USER  INTERFACE  (10/10) 
LEARNABILITY  (10/10) 
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BENCHMARK  CASES  DEVELOPED 
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CODE  LIMITATIONS 
LIMITED  LEARNING  VALUE 
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BENCHMARK  CASES  - SET  2 (CONT'D) 
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SCORING  GUIDELINES  DEVELOPED 
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GRIDGEN:  TFI  AND  EXTENSIVE  OPTIONS 
(SCORE  = 7+2  = 9 PTS) 


CANDIDATE  CODES  IDENTIFIED 
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CODE  EVALUATION 
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• DOCUMENTATION  AND  TRAINING  AVAILABLE 


1080 


• NASA  MSFC  EXPERIENCE  WITH  GENIE 

• ROCKETDYNE  EXPERIENCE  WITH  PATRAN  AND  IGB 

AVAILABILITY  DURING  PLANNED  EVALUATION  TIMEFRAME 
AVAILABLE  VERSIONS  "FROZEN"  FOR  PHASE  2 EVALUATION 


PHASE  2 - PRELIMINARY  EVALUATION 
MATRIX  DEVELOPED 
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PHASE  2 - PRELIMINARY  EVALUATION  MATRIX 

(Page  2 of  2) 


r 


if) 

O 

3 

CC 

O LO 

inoo-»-*-omcMcoco 

7 

8 
10 

- 

10 

10 

10 

1442.0 

PATRAN 

cm 

00C0LriTtl^r^-*M-M“CM 

0)0  10 

O 

7 

7 

5 

1072.0 

GRIDGEfv 

o> 

coior-^oinncooj 

CO  O)  N 

CM 

10 

5 

5 

1395.0 

GENIE 

in  n 

'f—  in  cm 

- 

10 

7 

10 

1106.0 

IGB 

CM  h- 

(0(0  0 

- 

10 

7 

10 

1135.0 

EAGLE- 

VIEW 

00  O 

(MNlOlrtOl^tOi-OO 

lO  N N 

- 

10 

5 

5 

1036.0 

CDC 

^ o 

OOCDCDCMCOinuO'M-tn 

7 

0 

6 

CM 

1 

1 

5 

1180.0 

WEIGHT 

10.000 

5.000 

10.000 

10.000 

4.000 

4.000 

8.000 
10.000 

7.000 

8.000 
8.000 

o o o 
o o o 
o o o 
o o d 

5.000 

10.000 

6.000 

6.000 

CRITERIA 

GRID  CONTROL 

• CLUSTERING 

• LOCAL  CONTROL 

USABILITY 

• LEARNABILITY 
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• SAVE  INTERMEDIATE  STATE 

• GRID  DIAGNOSTICS 
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• SOURCE  CODE 

• RESPONSIVE  TO  NEEDS 
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COST  OF  SERVICES 
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PHASE  3 - FINAL  EVALUATION  COMPLETED 
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REAL  TIME  COLOR  CODED  GRID  DIAGNOSTICS 
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STRUCTURED  ADAPTIVE  GRID  GENERATION  / 9 9 5"  / //,  99  9 
USING  ALGEBRAIC  METHODS 

by 

Jiann-CherDg  Yang,  Bharat  K.  Soni 
NSF/ERCCFS 

Mississippi  State  University,  MS39762 
and 

R.  P.  Roger,  Stephen  C.  Chan 
Teledyne  Brown  Engineering,  Huntsville,  AL  35807 

ABSTRACT 

The  accuracy  of  the  numerical  algorithm  depends  not  only  on  the  formal  order  of  approximation 
but  also  on  the  distribution  of  grid  points  in  the  computational  domain.  Grid  adaptation  is  a proce- 
dure which  allows  optimal  grid  redistribution  as  the  solution  progresses.  It  offers  the  prospect  of 
accurate  flow  field  simulations  without  the  use  of  execessively  time,  computationally  expensive, 
grid. 

Grid  adaptive  schemes  are  divided  into  two  basic  categories:  differentiall  and  algebraic.  The  dif- 
ferential method  is  based  on  a variational  approach  where  a function  which  contains  a measure  of 
grid  smoothness,  orthogonality  and  volume  variation  is  minimized  by  using  a variational  principle. 
This  approach  provideds  a solid  mathematical  basis  for  the  adaptive  method,  but  the  Euler-La- 
grange  equations  must  be  solved  in  addition  to  the  original  governing  equations.  On  the  other  hand, 
the  algebraic  method  requires  much  less  computational  effort,  but  the  grid  may  not  be  smooth.  The 
algebraic  techniques  are  based  on  devising  an  algorithm  where  the  grid  movement  is  governed  by 
estimates  of  the  local  error  in  the  numerical  solution.  This  is  achieved  by  requiring  the  points  in  the 
large  error  regions  to  attract  other  points  and  points  in  the  low  error  region  to  repel  other  points. 

The  development  of  a fast,  efficient,  and  robust  algebraic  adaptive  algorithm  for  structured  flow 
simulation  applications  is  presented.  This  development  is  accomplished  in  a three  step  process.  The 
first  step  is  to  define  an  adaptive  weighting  mesh  (distribution  mesh)  [1][2]  on  the  basis  of  the  equi- 
distribution  law  applied  to  the  flow  field  solution.  The  second,  and  probably  the  most  crucial  step, 
is  to  redistribute  grid  points  in  the  computational  domain  according  to  the  aforementioned  weighting 
mesh.  The  third  and  the  last  step  is  to  reevalute  the  flow  property  by  an  appropriate  search/interpolate 
scheme  at  the  new  grid  locations. 

The  adaptive  weighting  mesh  provides  the  information  on  the  desired  concentration  of  points 
to  the  grid  redistribution  scheme.  The  evaluation  of  the  weighting  mesh  is  accomplished  by  utilizing 
the  weight  function  representing  the  solution  variation  and  the  equidistribution  law  [3].  The  selec- 
tion of  the  weight  function  plays  a key  role  in  grid  adaptation  [4]  A new  weight  function  utilizing 
a properly  weighted  boolean  sum  of  various  flowfield  characteristics  is  defined.  The  redistribution 
scheme  is  developed  utilizing  Non-Uniform  Rational  B-Splines  (NURBS)  representation  [5],  The 
application  of  NURBS  representation  results  in  a well  distributed  smooth  grid  by  maintaining  the 
fidelity  of  the  geometry  associated  with  boundary  curves.  Several  algebraic  methods  are  applied  to 
smooth  and/or  nearly  orthogonalize  the  grid  lines.  Elliptic  solver  is  utilized  to  smooth  the  grid  lines 
if  there  are  grid  crossing. 

Various  computational  examples  of  practical  interest  are  presented  to  demonstrate  the  success 
of  these  methods.  (1)  Single  wedge  case.  (2)  Cylinder  case.  (3)  Airfoil  NACA  0012  case.  (4)  Viscous 
missile  case.  (5)  Slot  cooled  seeker  window  case.  (6)  Moving  shock  in  converge  inlet  case.  (7)  3-D 
single  wedge  case.  (8)  3-D  viscous  cone  case.  (9)  3-D  f 15  body  and  wing  case.  (10)  3-D  helicopter 
propeller  shock  prediction  case.  In  case  (l),(2),(4),(5),(6),(7),(8),and(9),  PARC2D  and  PARC3D  [6] 
are  used  as  the  flow  simulation  code.  In  case  (3),  UBI  [7]  is  used  as  the  flow  simulation  code,  and 
compare  to  Peter  M.  Goorjian’s  solution  [8].  In  case  (10),  MSUTC  [9]  is  used  as  the  flow  simulation 
code. 
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Figure  1. 

Use  NURBS  surface  to  generate  smooth  grid. 


Single  Wedge  Case 
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Figure  3. 

Single  wedge  grid  adaptation  process. 
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Figure  5. 

Use  B-Spline  curve  to  keep  sharp  corners. 


Cylinder  Case 
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Figure  6. 

Flow  around  cylinder  initial/adaptive  cases. 
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Figure  8. 

Airfoil  NACA0012  case  non-adaptive/adaptive  grid/solution, 
Mach  = 0.8,  a = 1.25. 


Viscous  missile 


Figure  11. 

Viscous  missile  case,  bow  shock  in  front  of  leading  edge. 
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Figure  13. 

Slot  cooled  seeker  window  case,  velocity  contour  near  interceptor. 


1115 


ifiWOTnd 


Figure  14. 

Adaptive  grid  for  the  moving  shock  in  the  converge  inlet. 
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GEOMETRY  I GEOMETRY 

30  initial  arid  1 30  second  adaptive  grid  in  1 direction 

for  single  Hedge  at  k*G  surface  ;1  for  single  Hedge  at  k-6  surface 


Figure  15. 

3-D  Single  wedge,  comparison  of  initial  and  adaptive  grid  & solution 


3-D  Viscous  Cone  Case 
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Figure  16. 

3-D  viscous  cone  case,  bow  shock  in  front  of  cone  head. 


3-D  F-15  Body  and  Wing 
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Figure  17. 

3-D  F-15  body  and  wing  case,  symmetric  surface. 
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3-D  initial 
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3-D  adaptive  grid 
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Figure  20. 

3-D  adaptive  grid  surfaces  for  the  helicopter  propeller  intersection. 
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Abstract 

The  objective  of  this  research  is  to  develop  an  efficient,  time-accurate  numerical 
algorithm  to  discretize  the  Navier-Stokes  equations  for  the  predictions  of  internal  one-, 
two-dimensional  and  axisymmetric  flows.  A generic,  efficient,  elliptic  adaptive  grid 
generator  is  implicitly  coupled  with  the  Lower-Upper  factorization  scheme  in  the 
development  of  ALUNS  computer  code.  The  calculations  of  one-dimensional  shock  tube 
wave  propagation  and  two-dimensional  shock  wave  capture,  wave-wave  interactions, 
shock  wave-boundary  interactions  show  that  the  developed  scheme  is  stable,  accurate  and 
extremely  robust.  The  adaptive  grid  generator  produced  a very  favorable  grid  network  by 
a grid  speed  technique.  This  generic  adaptive  grid  generator  is  also  applied  in  the  PARC 
and  FDNS  codes  and  the  computational  results  for  solid  rocket  nozzle  flowfield  and 
crystal  growth  modeling  by  those  codes  will  be  presented  in  the  conference,  too.  This 
research  work  is  being  supported  by  NASA/MSFC. 
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Shock  tube  pressure  solution 


1139 


Ill  I 111 


1140 


Pressure  contour  plot  inside  a two-dimensional  duct 
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Adapted  grid  of  a two-dimensional  duct 
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block  interfaces.  It  has  been  applied  to  generate  grids  for  various  complicated  tur- 
bomachinery geometries  (Figure  2),  as  well  as  to  rocket  and  missile  configurations 
(Figure  3). 
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o User  interaction  for  grid  manipulation 


X 

o 

cd 

o 

u 


z 

w 

o 

si 

•ph 

a? 


s 

3 


*h 

Sd 

13 

© 

© 

X 

*H 

«s 

l> 

t> 


d § 

1 I 

8 I 

A fa 


ft 

ft 

cd 

'd 

0 

ww 

*H 

pc 

• pH 

<1 

ft 

JO 

IS 

a 

•pH 

& 

GO 

s 

> 

© 

*d 

>> 

0 

X 

d 

tf 

0 

ft 

g 

© 

1 % 

X 

1 a 

cd 

a 

• fH 

S 

•F* 

fee 

hH 

•c 

© 

5 

O 

tf 

0 

g a 

| .2 

I -2 
I g. 

sj 
CD 

•PH 

1 


• fH 


a 

o 

•PH 

-M 

p2 

o 

& 

© 

'M 

a 

TJ 

•ph 

& 


a 

ja 

'o 

> 

© 

© 

_g 

«4H 

*H 

CD 


o 

PS 

u 

a 

hH 

ft 

H 

© 


I 


a 

o 

•pH 

cd 

© 

d 

© 

be 


© »ph 


So 

© 

a 

'o 

> 

© 

© 

*cd 


CD 

»H 

<2 

a 

© 

+* 

CD 

CD 

© 

• pH 


w fa 

•pN  *fH 


o 0 0 0 0 0 0 


1153 


April  21, 1993  CFD  Workshop 


TIGER 


Converts  various  alien  formats  for  geometry  definition 
2D  LE/TE  circle  fitting 


Applications 


PA 

GO 

3 § 
1 1 
£ i 

.a  § 

« o 

| JS 

>r  '8 

1 .2 

2 s 


(O  CO  09 

fn  Jm  5h 

P P P 

0 o o 

^ X 4= 
<m  c^q  <m 

1 l i 


H § 

0 £ 


c TT 

S| 

o ^ 

rH_  rH^ 

r>  ^ 

» p 

CO  TS 


-d  £ 
P n 

P S 

;£  w 

CO  Q 

3 O 


•w4  ^2 

a > 

P p 

n * 


^ o 
p p 

6 £ 

<N  K 


C Ss 

^ u 

^ -2 

g ?a 

H il  M 
0>  0 -3 
+2  -g  P 
d ,2  © 

-h 


fifif 

111 
CJ  O 9 
S 3 * 

cog 

S ’S  'Q 
3§a 

d 9 w 

8 ? >1 

| ◄ | 

i ss 

P 'M  63 

*a  h 

££  s 


2 £ 

CU  eg 


2 H 


April  21, 1993  CFD  Workshop 


O 0 0 0 


1156 


Rotor  Tip  Cutter  (Spherical) 
3D  Blade  LE/TE  Circle  Fitting 
CAD  Interface  (CAGI) 


Figure  1.  Tiger  Graphical  User  Interface 
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Figure  2.  Engine  Mockup 
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Towards  A Generalized  Computational  Fluid  Dynamics 

Technique  for  all  Mach  Numbers  Cq  ^ w- 


R.  W.  Walters,  D.  C.  Slack  and  A.  G.  Godfrey 


Statement  of  the  Problem  /79^>  / / 9 600 

Currently  there  exists  no  single  unified  approach  for  efficiently  and  accurately  solving 
computational  fluid  dynamics  (CFD)  problems  across  the  Mach  number  regime,  from  truly 
low  speed  incompressible  flows  to  hypersonic  speeds.  There  are  several  CFD  codes  that  have 
evolved  into  sophisticated  prediction  tools  with  a wide  variety  of  features  including  multi- 
block capabilities,  generalized  chemistry  and  thermodynamics  models  among  other  features. 
However,  as  these  codes  evolve,  the  demand  placed  on  the  end  user  also  increases  simply 
because  of  the  myriad  of  features  that  are  incorporated  into  these  codes.  In  order  for  a user  to 
be  able  to  solve  a wide  range  of  problems,  several  codes  may  be  needed  requiring  the  user  to 
be  familiar  with  the  intricacies  of  each  code  and  their  rather  complicated  input  files. 
Moreover,  the  cost  of  training  users  and  maintaining  several  codes  becomes  prohibitive. 


Objective  of  the  Work 

The  objective  of  the  current  work  is  to  extend  the  compressible,  characteristic-based, 
thermochemical  nonequilibrium  Navier-Stokes  code  GASP  to  very  low  speed  flows  and 
simultaneously  improve  convergence  at  all  speeds.  Before  this  work  began,  the  practical 
speed  range  of  GASP  was  Mach  numbers  on  the  order  of  0.1  and  higher.  In  addition,  a 
number  of  new  techniques  have  been  developed  for  more  accurate  physical  and  numerical 
modeling. 


Approach  Used 

The  primary  focus  has  been  on  the  development  of  optimal  preconditioning  techniques  for  the 
Euler  and  the  Navier-Stokes  equations  with  general  finite-rate  chemistry  models  and  both 
equilibrium  said  nonequilibrium  thermodynamics  models.  We  began  with  the  work  of  Van 
Leer,  Lee,  and  Roe  for  inviscid,  one-dimensional  perfect  gases  and  extended  their  approach  to 
include  three-dimensional  reacting  flows.  The  basic  steps  required  to  accomplish  this  task 
were  a transformation  to  steam-aligned  coordinates,  the  formulation  of  the  preconditioning 
matrix,  incorporation  into  both  explicit  and  implicit  temporal  integration  schemes,  and 
modification  of  the  numerical  flux  formulae.  In  addition,  we  improved  the  convergence  rate  of 
the  implicit  time  integration  schemes  in  GASP  through  the  use  of  inner  iteration  strategies 
and  thp  use  of  the  GMRES  (General  Minimized  RESisual)  which  belongs  to  the  class  of 
algorithms  referred  to  as  Krylov  subspace  iteration.  Finally,  we  significantly  improved  the 
practical  utility  of  GASP  through  the  addition  of  mesh  sequencing,  a technique  in  which 
computations  begin  on  a coarse  grid  and  get  interpolated  onto  successively  finer  grids. 


Conclusions  Relevant  to  Rocket  Propulsion 


The  fluid  dynamic  problems  of  interest  to  the  propulsion  community  involve  complex  flow 
physics  spanning  different  velocity  regimes  and  possibly  involving  chemical  reactions.  This 
/♦lagg  of  problems  results  in  widely  disparate  time  scales  causing  numerical  stiffness.  Even  in 
the  absence  of  chemical  reactions,  eigenvalue  stiffness  manifests  itself  at  transonic  and  very 
low  speed  flows  which  can  be  quantified  by  the  large  condition  number  of  the  system  and 
evidenced  by  slow  convergence  rates.  This  results  in  the  need  for  thorough  numerical 
analysis  and  subsequent  implementation  of  sophisticated  numerical  techniques  for  these 
difficult  yet  practical  problems.  As  a result  of  this  work,  we  have  been  able  to  extend  the 
range  of  applicablity  of  compressible  codes  to  very  low  speed  inviscid  flows  (M=.001)  and 
reacting  flows.  Our  work  now  centers  on  the  extension  to  viscous  flows. 
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Table  1:  AeroSoft  Phase  II  Timetable 
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PRECONDITIONING  FOR  THE  NAVIER-STOKES  EQUATIONS 
WITH  FINITE-RATE  CHEMISTRY 


Andrew  G.  Godfrey 

Aero  Thermo  Technology 
6703  Odyssey  Drive,  Suite  303 
Huntsville,  Alabama  35806 


f. 


ABSTRACT 


The  extension  of  Van  Leer’s  preconditioning  procedure  to  generalized  finite-rate  chem- 
istry is  discussed.  Application  to  viscous  flow  is  begun  with  the  proper  preconditioning 
matrix  for  the  one-dimensional  Navier-Stokes  equations.  Eigenvalue  stiffness  is  resolved 
and  convergence-rate  acceleration  is  demonstrated  over  the  entire  Mach-number  range  from 
nearly  stagnant  flow  to  hypersonic.  Specific  benefits  are  realized  at  the  low  and  transonic 
flow  speeds  typical  of  complete  propulsion-system  simulations.  The  extended  precondition- 
ing matrix  necessarily  accounts  for  both  thermal  and  chemical  non-equilibrium.  Numerical 
analysis  reveals  the  possible  theoretical  improvements  from  using  a preconditioner  for  all 
Mach  number  regimes.  Numerical  results  confirm  the  expectations  from  the  numerical 
analysis.  Representative  test  cases  include  flows  with  previously  troublesome  embedded 
high-condition-number  areas. 

Van  Leer,  Lee,  and  Roe  recently  developed  an  optimal,  analytic  preconditioning  tech- 
nique to  reduce  eigenvalue  stiffness  over  the  full  Mach-number  range.  By  multiplying 
the  flux-balance  residual  with  the  preconditioning  matrix,  the  acoustic  wave  speeds  are 
scaled  so  that  all  waves  propagate  at  the  same  rate,  an  essential  property  to  eliminate 
inherent  eigenvalue  stiffness.  This  session  discusses  a synthesis  of  the  thermo-chemical 
non-equilibrium  flux-splitting  developed  by  Grossman  and  Cinnella  and  the  characteris- 
tic wave  preconditioning  of  Van  Leer  into  a powerful  tool  for  implicitly  solving  two  and 
three-dimensional  flows  with  generalized  finite-rate  chemistry. 

For  finite-rate  chemistry,  the  state  vector  of  unknowns  is  variable  in  length.  Therefore, 
the  preconditioning  matrix  extended  to  generalized  finite-rate  chemistry  must  accommo- 
date a flexible  system  of  moving  waves.  Fortunately,  no  new  kind  of  wave  appears  in  the 
system.  The  only  existing  waves  are  entropy  and  vorticity  waves,  which  move  with  the 
fluid,  and  acoustic  waves,  which  propagate  in  Mach-number  dependent  directions.  The 
non-equilibrium  vibrational  energies  and  species  densities  in  the  unknown  state  vector  act 
strictly  as  convective  waves.  The  essential  concept  for  extending  the  preconditioning  to 
generalized  chemistry  models  is  determining  the  differential  variables  which  symmetrize 
the  flux  Jacobians.  The  extension  is  then  straight-forward. 

This  algorithm  research  effort  will  be  released  in  a future  version  of  the  production- 
level  computational  code  coined  the  General  Aerodynamic  Simulation  Program  (GASP), 
developed  by  Walters,  Slack  and  McGrory. 
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H.M.  Shang  and  Y.S.  Chen 
Engineering  Sciences,  Inc.,  Huntsville,  AL 


Z.J.  Chen  and  C.P.  Chen 
University  of  Alabama  in  Huntsville,  AL 


T.S.  Wang 

NASA  MSFC,  Huntsville,  AL 


ABSTRACT 


Combustion  processes  in  rocket  propulsion  systems  are  characterized  by 
the  existence  of  multiple,  vastly  differing  time  and  length  scales,  as  well  as 
flow-speeds  at  wide  variation  of  Mach  numbers.  The  chemical  kinetics 
processes  in  the  highly  active  reaction  zone  are  characterized  by  much  more 
smaller  scales  compared  to  fluid  convective  and  diffusive  time  scales.  An 
operator  splitting  procedure  for  transient  finite  rate  chemistry  problems  has 
been  developed  using  a pressure  based  method,  which  can  be  applied  to  all 
speed  flows  without  difficulties.  The  splitting  of  chemical  kinetics  terms 
formed  the  fluid-mechanical  terms  of  the  species  equation  ameliorated  the 
difficulties  associated  with  the  disparate  time  scales  and  stiffness  in  the  set 
of  equations  which  describes  highly  exothermic  combustion.  A combined 
efficient  ordinary  differential  equations  (O.D.E.)  solver  was  used  to  integrate 
the  effective  chemical  source  terms  over  the  residence  time  at  each  grid  cell. 
One  and  two  dimensional  reacting  flow  situations  were  carried  out  to 
demonstrate  and  verify  the  current  procedure.  Different  chemical  kinetics 
with  different  degrees  of  nonlinearity  have  also  been  incorporated  to  test  the 
robustness  and  generality  of  the  proposed  method. 
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Three-Dimensional  Navier-Stokes  Analysis  and  Redesign  of  an 
Imbedded  Bellmouth  Nozzle  in  a Turbine  Cascade  Inlet  Section 


P.  W.  Giel,  Sverdrup  Technology,  Inc.  (216)  826-6686 

J.  R.  Sirbaugh  NASA  LeRC  Group 

2001  Aerospace  Parkway 
Brook  Park,  OH  44142 


Verification  of  proposed  turbopump  blading  performance  will  involve  evaluation 
of  candidate  blades  in  cascade  test  facilities.  It  is  necessary  to  be  able  to  predict  the  flow 
Fields  within  these  cascades  for  the  results  to  be  applicable  to  actual  engine 
environments.  This  work  presents  the  results  of  a study  to  predict  the  flow  field  for  the 
NASA  Lewis  Transonic  Turbine  Blade  Cascade  Facility,  which  is  similar  to  those  used 
to  evaluate  rocket  propulsion  turbines.  A pitchwise  non-uniform  total  pressure 
distribution  was  observed  at  the  blade  row  leading  edge  plane.  A CFD  analysis  was 
used  to  show  that  the  cause  of  the  flow  non-uniformity  was  a pair  of  vortices  that 
originated  in  an  embedded  bellmouth  inlet.  Further  CFD  analysis  was  used  to  verify 
that  a redesigned  inlet  section  resulted  in  a flow  with  acceptable  uniformity. 

A computational  analysis  was  chosen  because  physical  accessibility  to  the  inlet 
section  was  limited,  and  because  a computational  approach  also  allows  one  to  examine 
design  changes  cheaper  and  more  quickly  than  an  experimental  approach  would.  The 
PARC  code,  a general  purpose,  three-dimensional,  Navier-Stokes  code  with  multi-block 
solution  capability,  was  chosen  for  the  present  study.  Results  are  presented  detailing  the 
computational  requirements  needed  to  accurately  predict  flows  of  this  nature. 

Calculations  of  the  original  geometry  showed  total  pressure  loss  regions  consistent 
in  strength  and  in  location  to  experimental  measurements.  An  examination  of  the 
results  shows  that  the  distortions  are  caused  by  a pair  of  vortices  that  originate  as  a 
result  of  the  interaction  of  the  flow  with  the  imbedded  bellmouth.  Computations  were 
performed  for  an  inlet  geometry  which  eliminated  the  imbedded  bellmouth  by  bridging 
the  region  between  it  and  the  upstream  wall.  This  analysis  indicated  that  eliminating 
the  imbedded  bellmouth  eliminates  the  troublesome  pair  of  vortices,  resulting  in  a flow 
with  much  greater  pitchwise  uniformity. 
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Transonic  Turbine  Blade  Cascade  Inlet  Analysis 
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off-design  incidence  case 


Transonic  Turbine  Blade  Cascade  Inlet  Analysis  - 
Overall  View  of  Cascade  Inlet  Section 
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Transonic  Turbine  Blade  Cascade  Inlet  Analysis  - 
Preliminary  Experimental  Measurements 
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Transonic  Turbine  Blade  Cascade  Inlet  Analysis  - 
Code  Descriptions 
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Second-order  accurate  finite-differences 
Blended  second-  and  fourth-difference  artificial  dissipation 
Four-stage  Runge-Kutta  solution  algorithm 
Baldwin-Lomax  or  Cebeci-Smith  algebraic  turbulence  models 


Transonic  Turbine  Blade  Cascade  Inlet  Analysis  - 
Original  Geometry  and  Grids 
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Transonic  Turbine  Blade  Cascade  Inlet  Analysis  - 
Calculations  of  Original  Geometry  (uniform  exit) 


Transonic  Turbine  Blade  Cascade  Inlet  Analysis  - 
RVC3D  Isolated  Blade  Calculations  - Loading 
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Transonic  Turbine  Blade  Cascade  Inlet  Analysis  - 
Non-uniform  Exit  Pressure  Boundary  Condition 


Transonic  Turbine  Blade  Cascade  Inlet  Analysis  - 
Calculations  of  Original  Geometry  (non-uniform  exit) 


Transonic  Turbine  Blade  Cascade  Inlet  Analysis  - 
Calculations  of  Original  Geometry  (non-uniform  exit) 


Transonic  Turbine  Blade  Cascade  Inlet  Analysis  - 
Calculations  of  Original  Geometry  (non-uniform  exit) 


Transonic  Turbine  Blade  Cascade  Inlet  Analysis  - 
Modified  Geometry  and  Grids 
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Transonic  Turbine  Blade  Cascade  Inlet  Analysis 
Upstream  8:1  Contraction;  geometry  and  grid 
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Transonic  Turbine  Blade  Cascade  Inlet  Analysis  - 
Upstream  8:1  Contraction;  preliminary  results 
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Transonic  Turbine  Blade  Cascade  Inlet  Analysis  - 

Calculations  of  Modified  Geometry  with 

Unstream  8:1  Contraction  (non-uniform  exit) 
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Abstract 

Prediction  of  Incidence  and  Surface  Roughness 

Turbine  Performance 

R.  J.  Boyle 

NASA  Lewis  Research  Center 
Cleveland,  OH  44135 

The  results  of  a Navier-Stokes  analysis  for  predicting  the 
change  in  turbine  efficiency  due  to  a change  in  either 
incidence  or  surface  roughness  is  discussed.  It  was 
experimentally  determined  by  Boynton,  Tabibzadeh,  and  Hudson 
that  polishing  the  SSME  high  pressure  fuel  turbine  blades 
improved  turbine  efficiency  by  about  2 points  over  a wide 
range  of  operating  conditions.  These  conditions  encompassed 
the  range  of  incidence  seen  by  the  turbine  blading  during 
flight.  It  is  also  necessary  to  be  able  to  predict  turbine 
performance  at  various  operating  points  for  future  rocket 
turbopump  applications.  The  code  RVCQ3D,  developed  by  Rod 
Chima,  was  used  to  determine  the  effects  of  changes  in 
incidence  angle  on  turbine  blade  row  efficiency.  The  midspan 
Navier-Stokes  results  were  used  in  conjunction  with  an 
inviscid  flow  analysis  code  to  predict  the  efficiency  of  the 
two  stage  SSME  over  a wide  range  of  operating  conditions  for 
smooth  and  rough  turbine  blades.  The  use  of  the  Navier- 
Stokes  analysis  to  predict  changes  in  turbine  efficiency  due 
to  variation  in  incidence  angles  was  found  to  be  superior  to 
other  incidence  loss  correlations  available  in  the 
literature.  The  sensitivity  of  the  Navier-Stokes  results  to 
grid  parameters  is  discussed. 

The  effects  of  the  surface  roughness  were  accounted  for 
using  the  Cebeci-Chang  rough  wall  turbulence  model.  This 
model  was  implimented  in  the  code  RVCQ3D.  The  implimentation 
of  this  model  for  predicting  the  change  in  efficiency  is 
also  discussed. 
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INCIDENCE  ar>d  SURFACE  ROUGHNESS 
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COMPUTATIONAL  APPROACH 
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VISCOUS  - RVCQ3D  - MIDSPAN 
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SECONDARY  FLOW 
INCIDENCE 


VISCOUS  SOLVER 
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MODIFIED  FOR  ROUGHNESS  EFFECTS 

CEBECI-CHANG  MODEL 
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MODEL  IN  PERFORMANCE  CALCULATION 
(MTSB) 


INCIDENCE  LOSS  COMPARISONS 
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40%  SPEED  LINE  vs  PRESSURE  RATIO 
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TRANSONIC  FLOW  CONDITIONS 


SENSITIVITY  OF  ROTOR  LOSS  TO  GRID  PARAMETERS 
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EFFICIENCY  FOR  SMOOTH  BLADED  SSME  HPFT 
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EFFICIENCY  FOR  SMOOTH  BLADEDED  SSME  HPFT 


NORMALIZED  WHEEL  SPEED,  Um/VinKAL 


LOSS  PREDICTIONS  - NAVIER-STOKES  & CORRELATIONS 

SSME  HPFT  FIRST  STAGE  ROTOR 


INCIDENCE  ANGLE,  DEG 


LOSS  PREDICTIONS  - NAVIER-STOKES  & CORRELATIONS 
SSME  HPFT  SECOND  STAGE  STATOR 
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INCIDENCE  ANGLE,  DEG. 


EFFECT  OF  OPTIMUM  INCIDENCE  ANGLE  ON  LOSS 

SSME  FIRST  STAGE  ROTOR 
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EFFECT  OF  OPTIMUM  INCIDENCE  ANGLE  ON  LOSS 
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INCIDENCE  ANGLE,  DEG. 


EFFICIENCY  OF  SMOOTH  BLADED  SSME  HPFT 
NAVIER-STOKES  DERIVED  INCIDENCE  LOSS  MODEL 
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INCIDENCE  LOSS  CONCLUSIONS 
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IMPROVED  INCIDENCE  LOSS  CORRELATION  VIA 
VARIATION  OF  NAVIER-STOKES  PARAMETERS 

NOT  ALWAYS  CONVENIENT  IN  EXPERIMENTS 
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NEAR  WALL  DAMPING  (A")  PART  OF  TURBULENT 
EDDY  VISCOSITY  MODEL 
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Abstract  of  a proposed  presentation  at  the  Workshop  / 

for  CFD  Applications  in  Rocket  Propulsion  to  be  held 
at  NASA  Marshall  Space  Flight  Center,  AL,  April  20-22,  1993.  ^ f 
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Three-dimensional  unsteady  flow  calculations  in  an 
advanced  Gas  Generator  turbine 


Akil  A.  Rangwalla 
MCAT  Institute,  Moffett  Field,  CA 

Abstract 


This  paper  deals  with  the  application  of  a three-dimensional,  unsteady  Navier-Stokes 
code  for  predicting  the  unsteady  flow  in  a single  stage  of  an  advanced  gas  generator  turbine. 
The  numerical  method  solves  the  three-dimensional  thin-layer  Navier-Stokes  equations, 
using  a system  of  overlaid  grids,  which  allow  for  relative  motion  between  the  rotor  and 
stator  airfoils.  Results  in  the  form  of  time  averaged  pressures  and  pressure  amplitudes  on 
the  airfoil  surfaces  will  be  shown.  In  addition,  instantaneous  contours  of  pressure,  mach 
number  etc.  will  be  presented  in  order  to  provide  a greater  understanding  of  the  inviscid 
as  well  as  the  viscous  aspects  of  the  flowfield.  Also,  relevant  secondary  flow  features  such 
as  cross-plane  velocity  vectors  and  total  pressure  contours  will  be  presented.  Prior  work  in 
two- dimensions  has  indicated  that  for  the  advanced  designs,  the  unsteady  interactions  can 
play  a significant  role  in  turbine  performance.  These  interactions  affect  not  only  the  stage 
efficiency  but  can  substantially  alter  the  time-averaged  features  of  the  flow.  This  work  is 
a natural  extension  of  the  work  done  in  two- dimensions  and  hopes  to  address  some  of  the 
issues  raised  by  the  two-dimensional  calculations.  These  calculations  are  being  performed 
as  an  integral  part  of  an  actual  design  process  and  demonstrate  the  value  of  unsteady 
rotor-stator  interaction  calculations  in  the  design  of  turbomachines. 
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THREE-DIMENSIONAL  UNSTEADY  FLOW  CALCULATIONS 
FOR  AN  ADVANCED  GAS  GENERATOR  TURBINE 

(Preliminary  results) 
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OXIDIZER  TURBINE  BASELINE  DESIGN 
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Background  contd. .. 
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THEIR  USE  IN  THE  DESIGN  PROCESS  OF  A TURBOMACHINE 
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Turbulent  eddy  viscosity  computed  using  Baldwin-Lomax 


Boundary  conditions 
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Study  of  accuracy 
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Nonlinear  rotor- stator  interactions 


Figure  Hoi -streak  calculation:  original  and  adapted  grids  for  the  stator 
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Figure  Hot  streak  calculation:  original  and  adapted  grids  for  the  rotor 
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MACH  CONTOURS  OBTAINED  FROM  A STRUCTURED  METHOD  (STAGE-2) 


oooooooooaaootr. 


1301 


ORIGINAL  page  IS 

OF  POOR  QU&OT 


Preliminary  3D  Results  for  the  GGOT 
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Geometry  Rescaling 
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Number  of  stator  blades 
Number  of  rotor  blades 


SURFACE  GRID  FOR  THE  GGOT 
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Turbine  operating  conditions 
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INSTANTANEOUS  MACH  NUMBERS  IN  THE  GGOT 
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Pressure  Coefficient  (Pstatic/Ptotal) 


1312 


1.4  -1.2  -1.0  -0.8  -0.6  -0.4 

Axial  Distance  (x) 


-0.2 


1313 


1314 


T/P 


1315 


MACH  CONTOURS  AT  20%  SPAN 
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MACH  CONTOURS  AT  MID-SPAN 
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NUMERICAL  SIMULATION  OF  STEADY  AND  UNSTEADY  VISCOUS  FLOW  IN 
TURBOMACHINERY  USING  PRESSURE  BASED  ALGORITHM 


B.  LAKSHMINARAYANA,  Y.  HO  and  A.  BASSON 


The  Pennsylvania  State  University 
Department  of  Aerospace  Engineering 
University  Park,  PA  16802 


The  objective  of  this  research  is  to  simulate  steady  and  unsteady  viscous  flows,  including 
rotor/stator  interaction  and  tip  clearance  effects  in  turbomachinery. 

The  numerical  formulation  for  steady  flow  developed  here  includes  an  efficient  grid 
generation  scheme,  particularly  suited  to  computational  grids  for  the  analysis  of  turbulent 
turbomachinery  flows  and  tip  clearance  flows,  and  a semi-implicit,  pressure-based  computational 
fluid  dynamics  scheme  that  directly  includes  artificial  dissipation,  and  is  applicable  to  both  viscous 
and  inviscid  flows.  The  values  of  these  artificial  dissipation  is  optimized  to  achieve  accuracy  and 
convergency  in  the  solution.  The  numerical  model  is  used  to  investigate  the  structure  of  tip  clearance 
flows  in  a turbine  nozzle.  The  structure  of  leakage  flow  is  captured  accurately,  including  blade-to- 
blade  variation  of  all  three  velocity  components,  pitch  and  yaw  angles,  losses  and  blade  static 
pressures  in  the  tip  clearance  region.  The  simulation  also  includes  evaluation  of  such  quantities  of 
leakage  mass  flow,  vortex  strength,  losses,  dominant  leakage  flow  regions  and  the  spanwise  extent 
affected  by  the  leakage  flow.  It  is  demonstrated,  through  optimization  of  grid  size  and  artificial 
dissipation,  that  the  tip  clearance  flow  field  can  be  captured  accurately. 

The  above  numerical  formulation  was  modified  to  incorporate  time  accurate  solutions.  An 
inner  loop  iteration  scheme  is  used  at  each  time  step  to  account  for  the  non-linear  effects.  The 
computation  of  unsteady  flow  through  a flat  plate  cascade  subjected  to  a transverse  gust  revels  that 
the  choice  of  grid  spacing  and  the  amount  of  artificial  dissipation  is  critical  for  accurate  prediction 
of  unsteady  phenomena.  The  rotor-stator  interaction  problem  is  simulated  by  starting  the 
computation  upstream  of  the  stator,  and  the  upstream  rotor  wake  is  specified  from  the  experimental 
data.  The  results  show  that  the  stator  potential  effects  have  appreciable  influence  on  the  upstream 
rotor  wake.  The  predicted  unsteady  wake  profiles  are  compared  with  the  available  experimental  data 
and  the  agreement  is  good,  the  numerical  results  are  interpreted  to  draw  conclusions  on  the  unsteady 
wake  transport  mechanism  in  the  blade  passage. 
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NUMERICAL  METHOD 
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blade  compared  to  the  flat  plate  cascade 
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Tip  clearance  flow  through  a turbine  cascade 
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Reduced  Frequency  (cjC/2U) 
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Fraction  of  Axial  Chord 
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UNSTEADY  VISCOUS  OVER  AN  AIRFOIL 

The  unsteadiness  is  generated  by  two  flapping  foils  oscillating  about  one 
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UNSTEADY  VISCOUS  OVER  AN  AIRFOIL 
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UNSTEADY  VISCOUS  OVER  AN  AIRFOIL 
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X/C=0.61 1 X/C=0.784  X/C=0.972 

PRESSURE  SIDE 


AMPLITUDE  OF  VELOCITY  FOR  HARMONIC  n= 


PHASE  OF  VELOCITY  FOR  HARMONIC  n= 


S3E 


TIP  CLEARANCE  MODELLING 
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3D  STEADY  FLOW  IN  THE  TIP  CLEARANCE  REGION 

FLOW  OF  A TURBINE 
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3D  STEADY  FLOW  IN  THE  TIP  CLEARANCE  REGION 

FLOW  OF  A TURBINE 
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good  tip  gap  resolution 


EMBEDDED  H-GRID  TOPOLOGY 


Grid  Inside  Gap 


COMPARISON  OF  MEASURED  AND  PREDICTED 
PITCHWISE  FLOW  ANGLES 

Axial  Position  (%  Axial  Chord  from  Leading  Edge)  : 60. 
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Numerical 
Bindon  (1990) 


COMPARISON  OF  MEASURED  AND  PREDICTED 
PITCHWISE  FLOW  ANGLES 


Pitchwise  Flow  Angle  (deg) 
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Numerical 
Bindon  (1990) 
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Numerical 
Bindon  (1990) 


COMPARISON  OF  MEASURED  AND  COMPUTED  SPANWISE  FLO 
ANGLES  (POSITIVE  ANGLE  WHEN  FLOW  DIRECTED  TOWARD 
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Bindon  (1990) 


CONTOURS  OF  LOSS  COEFFICIENT  (cj 
PREDICTED  AND  MEASURED 
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Pitch  wise  Pitch  wise 


DEVELOPMENT  OF  MASS  AVERAGED  LOSS 
COEFFICIENT  (ct) 
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Axial  Chord 


omm'iM  page  m 

OF  POOR  QUA14TY 


(c)  Middle  of  Tip  Gap  (computed)\\l  (d)  ^ ^ Tjp  GaD 


VELOCITY  VECTORS  AT  CONSTANT  SPANWISE 
POSITIONS  PROJECTED  ON  TO  A BLADE-TO-BLADE 

(XY)  PLANE 
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VELOCITY  VECTORS  AT  CONSTANT  SPANWISE 
POSmONS  PROJECTED  ON  TO  A BLADE-TO-BLADE 
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VELOCITY  VECTORS  AT  CONSTANT  SPANWISE 
POSITIONS  PROJECTED  ON  TO  A BLADE-TO-BLADE 

(XY)  PLANE 
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Midspan  Exit  Velocity 


Secondary  Vorticity  Vectors 
Z=  0.50 


COMPUTED  SECONDARY  VORTICITY  VECTORS 

[ ( ] 
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CONCLUSIONS 

• Geometric  series  distribution  scheme 

Good  control  over  the  boundary  points 

• Algebraic-elliptic  grid  generation  scheme 

Good  control  over  clustering 
Good  control  over  orthogonality 
Enhanced  stability  of  elliptic  generation 

• Embedded  H-grid 

Single  block  discretization  for  tip  clearance  cases 
No  Modification  of  blade  tip  shape  required 
Retains  H-grid  connectivity  pattern 

• Effect  of  artificial  dissipation 

Numerical  accuracy  vs.  convergence  rate 
Minimum  artificial  dissipation  should  be  used 
• Modelling  of  tip  clearance  flows 

Major  physical  phenomena  captured 

Location  of  leakage  vortex,  pitchwise  and  spanwise  angles,  losses, 
static  pressures  predicted  accurately 
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velocity  profiles  and  blade  pressure  agree  with  the  measurements, 
including  the  magnitude  and  phase  angle. 
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GGOT  TOTAL  PRESSURE  LOSS  CONTROL  CONCEPT  EVALUATION^  ^ ^ ^ ' 

R.  F.  BLUMENTHAL 
AEROJET  PROPULSION  DIVISION 
SACRAMENTO,  CA 


Total  pressure  loss  is  one  of  the  most  important  parameters  in  the  design  of  a 
turbine.  This  parameter  effects  not  only  the  turbine  performance,  but  consequently  the 
engine  power  balance  and  engine  performance.  Computational  Fluid  Dynamics  (CFD) 
can  be  an  effective  tool  in  predicting  turbine  total  pressure  loss,  and  also  for  performing 
sensitivity  studies  to  achieve  an  optimal  design  with  respect  to  pressure  loss.  In  the 
present  study,  the  AEROVISC  code  was  used  to  predict  the  total  pressure  loss  in  the 
Turbine  Technology  Team  Gas  Generator  Oxidizer  Turbine  (GGOT). 

The  objectives  in  this  study  are  two-fold.  It  is  first  necessary  to  determine  an 
optimal  methodology  in  predicting  total  pressure  loss.  The  type  of  grid,  grid  density  and 
distribution  are  parameters  which  may  effect  the  loss  prediction.  Also,  the  effect  of 
using  a standard  K-e  turbulence  model  with  wall  functions  versus  a two-layer  turbulence 
model  needs  to  be  investigated.  The  use  of  grid  embedding  to  resolve  areas  with  high 
flow  gradients  needs  to  be  explored.  The  second  objective  of  the  study  is  to  apply  the 
optimal  methodology  toward  evaluating  different  tip  leakage  control  concepts. 

The  approach  taken  in  this  study  was  as  follows: 

1)  A nominal  baseline  case  was  run  (baseline  grid  with  standard  wall 

functions) 

a)  Grid  parametrics  were  performed  on  grid  density 

b)  Grid  embedding  was  applied  to  the  rotor  leading  and  trailing  edges, 
and  in  the  tip  region. 

c)  Evaluation  of  a two-layer  turbulence  model  (in  progress) 

Each  of  the  above  cases  were  assessed  in  terms  of  total  pressure  loss  in 
comparison  with  the  baseline  case,  and  in  terms  of  the  difference  in  secondary  flow 
resolution  in  comparison  with  the  baseline  case. 

2)  The  optimal  methodology  from  Step  1 is  applied  towards  evaluating 

different  tip  leakage  control  concepts  which  will  include 

a)  Hollow  rotor 

b)  Hollow  rotor  with  partitions  (labyrinth  seal  approach) 

c)  Hollow  rotor  with  partitions  and  suction-side  rotor  slots  (to  reduce 

fluid  impingement  angle) 

As  this  work  is  still  in  progress,  conclusions  are  not  available  at  this  time. 
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RBINE  BASEUNE  DESIGN 
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TURBULENT  WALL  (STATIONARY  W.R.T. 

RELATIVE  FRAME) 
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EXIT  STATIC  PRESSURE  USED  FOR  PRELIMINARY  PARAMETERICS 
BASED  ON  STATIC  PRESSURE  SPECIFIED  AT  ONE  NODE  IN  EXIT  PLANE 
(PS=215S  PSI  AT  CENTER  OF  EXIT  PLANE) 
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(PS=2155  PSI  AT  CENTER  OF  EXIT  PLANE) 
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STATUS:  MODEL  IS  CURRENTLY  RUNNING.  CONVERGENCE  DIFFICULTIES  HAVE  BEEN  ENCOUNTERED. 
PROBABLY  RELATED  TO  HIGH  ASPECT  RATIOS  ( - 25000  AT  BLADE  SURFACE  AT  MID-SPAN). 
SOLUTIONS  BEING  EXAMINED  INCLUDE  INCREASING  NODE  COUNT  IN  RADIAL  DIRECTION. 
RUNNING  CODE  IN  DOUBLE  PRECISION. 
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lATES  WILL  THEN  BE  ASSESSED  AT  THE  NOMINAL  CLEARANCE. 
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R=GAS  CONSTANT  FT/DEG  R 
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ABSTRACT 


The  Gas  Generator  Oxidizer  Turbine  (GGOT)  Blade  is  being  analyzed  by  various  investigators 
under  the  NASA  MSFC-sponsored  Turbine  Stage  Technology  Team  design  effort.  The  present  work 
concentrates  on  the  tip  clearance  region  flow  and  associated  losses;  however,  flow  details  for  the 
passage  region  are  also  obtained  in  the  simulations.  The  present  calculations  simulate  the  rotor 
blade  row  in  a rotating  reference  frame  with  the  appropriate  coriolis  and  centrifugal  acceleration  terms 
included  in  the  momentum  equations.  The  upstream  computational  boundary  is  located  about  one 
axial  chord  from  the  blade  leading  edge.  The  boundary  conditions  at  this  location  have  been 
determined  by  Pratt  & Whitney  using  an  Euler  analysis  without  the  vanes  to  obtain  approximately  the 
same  flow  profiles  at  the  rotor  as  were  obtained  with  the  Euler  stage  analysis  including  the  vanes. 
Inflow  boundary  layer  profiles  are  then  constructed  assuming  the  skin  friction  coefficient  at  both  the 
hub  and  the  casing.  The  downstream  computational  boundary  is  located  about  one  axial  chord  from 
the  blade  trailing  edge,  and  the  circumferentially  averaged  static  pressure  at  this  location  was  also 
obtained  from  the  P&W  Euler  analysis. 

Results  obtained  for  the  3-D  baseline  GGOT  geometry  at  the  full  scale  design  Reynolds  number 
show  a region  of  high  loss  in  the  region  near  the  casing.  Particle  traces  in  the  near  tip  region  show 
vortical  flow  behavior  of  the  fluid  which  passes  through  the  clearance  region  and  exits  at  the 
downstream  edge  of  the  gap.  In  an  effort  to  reduce  clearance  flow  losses,  the  mini-shroud  concept 
was  proposed  by  the  Pratt  & Whitney  design  team.  Calculations  were  performed  on  the  GGOT 
geometry  with  the  mini-shroud.  Results  of  these  calculations  indicate  that  the  mini-shroud  does  not 
significantly  affect  the  flow  in  the  passage  region,  and  although  the  tip  clearance  flow  is  different,  the 
mini-shroud  does  not  seem  to  prevent  the  above-mentioned  vortical  flow  behavior.  Since  both  flow 
distortion  and  total  pressure  losses  are  similar  for  both  geometries,  the  addition  of  the  mini-shroud 
does  not  seem  to  reduce  the  tip  clearance  flow  effects. 


This  work  was  supported  by  NASA  Marshall  Space  Flight  Center  under  Contract  NAS8-38S65. 
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APPROACH 


UTILIZE  SRA  MINT  CODE 
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• FULL  NAVIER-STOKES  EQUATIONS 

• NO-SLIP  WALL  BOUNDARY  EQUATIONS  WITH 
SUBLAYER  RESOLUTION 

• ALGEBRAIC  MIXING  LENGTH  TURBULENCE  MODEL 

• IMPLICIT  LINEARIZED  BLOCK  SOLVER  (ADI) 
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FLOW  PARAMETERS 


SUPPLIED  BY  P&W  DESIGN  TEAM 

CIRCUMFERENTIALLY  - AVERAGED  SPANWISE  DISTRIBUTIONS 
FROM  EULER  CODE 
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28  POINTS  FROM  HUB  TO  TIP 
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Baseline  GGOT  with  Clearance 
Blade  Surface  Pressure  (92.2%  span) 
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CONCLUSIONS 


• BOTH  GEOMETRIES  SHOW  SIMILAR  VORTICAL  FLOW  BEHAVIOR 

• TOTAL  PRESSURE  LOSSES  ARE  NOT  SIGNIFICANTLY 
DIFFERENT 

• ADDITION  OF  MINI-SHROUD  DOES  NOT  REDUCE  THE  TIP 
CLEARANCE  FLOW  EFFECTS 
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SUPERSONIC  FLOW  AND  SHOCK  FORMATION  IN  TURBINE  TIP  GAPS  ' 


Summary 


John  Moore 


Mechanical  Engineering  Department 
Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  Virginia  24061-0238 


P.  // 


Shock  tormation  due  to  overexpansion  of  supersonic  flow  at  the  inlet  to  the 
tip  clearance  gap  of  a turbomachine  has  been  studied. 

As  the  flow  enters  the  tip  gap,  it  accelerates  around  the  blade  pressure-side 
corner  creating  a region  of  minimum  static  pressure.  The  "free  streamline" 
separates  from  the  wall  at  the  corner;  and,  for  Mach  numbers  greater  than  about 
1.3,  it  curves  back  to  intersect  the  blade  tip.  At  this  point,  the  freestream  flow  is 
abruptly  turned  parallel  to  the  surface,  giving  rise  to  an  oblique  shock. 

The  results  are  consistent  with  compressible  sharp-edged  orifice  flow 
calculations  found  in  the  literature  and  with  the  theory  of  oblique  shock  wave 
formation  in  supersonic  flow  over  a wedge.  For  freestream  Mach  numbers  of  1 .4 
to  1.8,  wave  angles  are  43  to  54  degrees,  and  turning  angles  are  9 to  20  degrees; 
as  the  Mach  number  increases,  the  angle  of  turn  also  increases. 

It  appears  that  in  a turbine,  after  separating  from  the  inlet  corner,  the  flow 
reattaches  on  the  blade  tip  and  an  oblique  shock  is  formed  at  0.4-1. 4 tip  gap 
heights  into  the  clearance  gap.  The  resulting  shock-boundary  layer  interaction  may 
contribute  to  further  enhancement  of  already  high  heat  transfer  to  the  blade  tip  in 
this  region.  This  in  turn  could  lead  to  higher  blade  temperatures  and  adversely 
affect  blade  life  and  turbine  efficiency. 


Introduction 

Tip  leakage  flow  through  the  clearance  gaps  of  unshrouded  turbomachinery 
blades  is  known  to  cause  reductions  in  efficiency  and  performance  [Roelke,  1973, 
Hourmouziadis  and  Albrecht,  1987].  It  is  also  known  that  enhanced  heat  transfer 
to  the  tips  of  turbine  rotor  blades,  resulting  from  the  separation  and  reattachment 
of  the  leakage  flow,  can  be  a major  factor  in  determining  blade  life  in  high 
temperature  gas  turbines  [Moore  et  al.,  1989].  There  have  therefore  been  many 
recent  studies  of  flow  and  heat  transfer  in  tip  gaps  [Bindon,  1987,  Metzger  and 
Bunker,  1989,  Metzger,  Dunn,  and  Hah,  1991,  Yaras  et  al.,  1989].  But  most  of 
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these  studies  have  considered  only  incompressible  flow,  and  features  of 
compressible  flow  such  as  overexpansion  to  high  supersonic  Mach  numbers,  shock 
formation  within  the  tip  gap,  and  shock  boundary  layer  interaction  have  received 
little  attention.  It  is  important  that  this  compressible  flow  physics  be  understood 
sufficiently  weW  to  be  included  in  the  turbine  design  process. 

Consider,  for  example,  turbine  rotor  blades  in  gas  turbines  operating  with 
transonic  flow.  Around  the  airfoil,  the  flow  accelerates  to  supersonic  Mach 
numbers  near  the  suction  surface  in  regions  of  low  static  pressure.  Similarly  flow 
passing  through  the  tip  clearance  gap  will  accelerate  as  the  pressure  falls.  But  its 
path  is  not  as  smooth  as  that  around  the  airfoil  profile.  Efforts  are  made  to  reduce 
the  leakage  flow  by  using  sharp  corners,  for  example  on  the  pressure  side,  and 
cavities.  The  flow  may,  therefore,  overexpand  locally  and  exhibit  regions  of 
supersonic  flow  and  complex  compressible  flow  structure. 

In  an  attempt  to  shed  some  light  on  compressible  flow  development  in  tip 
gaps,  Henry  and  Moore  [Moore  et  al.,  1989]  made  a preliminary  study  using  a 
water  table  flow. simulation.  Using  the  hydraulic  analogy  between  free  surface 
liquid  flow  and  two-dimensional  compressible  flow,  they  gained  some  insight  into 
possible  Mach  number  distributions  and  shock  patterns  to  be  found  in  turbine  tip 
clearance  gaps.  It  appeared  that  overexpansion  of  flow  around  the  pressure  surface 
corner  leads  to  an  oblique  shock  wave  which  extends  from  the  inlet  corner  region 
to  the  shroud  wall  at  about  two  tip  gap  heights  into  the  tip  gap  and  then  reflects 
back  to  the  turbine  blade  tip.  For  tip  gap  exit  pressures  corresponding  to  blade 
suction  surface  Mach  numbers  greater  than  1 .0,  they  found  local  maximum  Mach 
numbers  within  the  tip  gap  in  the  range  of  about  1.5  to  1.8. 

Compressible  Orifice  Flows 

Although  little  was  found  in  the  literature  about  compressible  tip  gap  flows, 
much  research  has  been  done  on  compressible  orifice  flows.  Benson  and  Pool 
[1965]  describe  early  research  in  this  area. 

Benson  and  Pool  numerically  solved  the  equations  for  steady,  isentropic  flow 
of  air  through  a two-dimensional  slit.  They  then  compared  the  results  with 
Schlieren  photographs  and  interferograms.  Figure  1 shows  the  computed  free 
streamlines  for  various  back  pressures.  The  flows  range  from  incompressible  flow 
pb/p0  ~ 1.0,  M — 0.0,  to  sonic  flow,  pb/p0  = 0.5283,  M = 1.0,  to  choked  flow, 
Ph/p*  = 0.0389,  M = 2.77. 

Consider  flow  accelerating  along  the  orifice  wall  to  the  sharp  corner.  If  the 
back  pressure  is  low  enough  to  cause  supersonic  flow,  this  flow  will  undergo  a 
Prandtl-Meyer  expansion  at  the  corner  to  the  freestream  Mach  number.  As  a 
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result  the  flow  turns  through  the  corresponding  expansion  angle.  The  free 
streamline  then  continues  to  turn  through  a further  ninety  degrees  to  the  point  of 
its  maximum  (downward)  slope.  Subsequently  the  jet  reaches  a maximum  width 
followed  by  a contraction  (as  the  free  streamline  turns  upward  again). 

For  subsonic  and  sonic  flow,  the  jet  simply  contracts  to  an  asymptotic  jet 
width.  This  width  is  ir/(7r+2)  = 0.61  of  the  slit  width  for  incompressible  flow 
and  0.74  of  the  slit  width  for  sonic  flow,  as  seen  in  Fig.  1.  With  supersonic  flow, 
the  jet  width  contracts  and  then  expands.  At  a freestream  Mach  number  of  about 
1.3  (pb/po  ~ 0.36),  the  maximum  width  of  the  jet  is  equal  to  the  slit  width.  At 
choked  flow,  the  maximum  width  is  about  five  times  the  slit  width. 

Norwood  [1961]  performed  a computational  study  similar  to  that  of  Benson 
and  Pool.  He  concentrated  mostly  on  flow  near  the  slit.  He  also  performed 
experimental  work  on  jet  reattachment  in  two-dimensional  models  of  a flapper 
valve.  The  geometry  of  these  was  like  that  of  the  slit  in  Fig.  1 except  that 
downstream  of  the  minimum  area  a straight  wall  was  inclined  downwards  at  an 
angle  of  22-1/2  degrees  to  the  horizontal.  This  effectively  produced  a tip  gap  with 
a sharp  inlet  corner  and  a linearly  increasing  height.  Norwood  visualized  the  flow 
with  shadowgraph  pictures. 

At  low  upstream  stagnation  pressures  (pb/p0  > ~ 0.23),  the  jet  in 
Norwood’s  models  followed  the  horizontal  wall  like  the  subsonic  and  just  super- 
sonic flows  in  Figure  1 . Further  increase  in  the  upstream  pressure  (pb/p0  < ~ 
0.23),  however,  caused  the  jet  to  jump  to  the  inclined  wall.  Norwood  noted  that 
in  his  two-dimensional  models  the  flow  jumped  at  a relatively  constant  pressure 
ratio.  This  suggested  that  the  only  characteristic  parameter  for  the  flow 
development  was  the  Mach  number.  It  also  indicated  that  the  phenomenon  of 
reattachment  on  the  inclined  wall  was,  initially  at  least,  a compressibility  effect 
depending  on  the  Mach  number  rather  than  a frictional  effect. 

In  a shadowgraph  picture  of  flow  along  the  inclined  wall,  Norwood  observed 
two  features  of  interest  in  the  present  study.  The  first  was  "a  teardrop  shaped 
region,"  or  "bubble,"  at  the  edge  of  the  orifice  where  the  flow  is  rapidly  accelerat- 
ing and  the  streamlines  are  highly  curved.  Norwood  argued  that  the  pressure  in 
the  bubble  was  very  low  because  of  the  entrainment  of  the  air  inside  into  the  main 
stream.  The  second  interesting  feature  was  an  oblique  shock  just  downstream  of 
the  throat  starting  on  the  inclined  wall  at  the  end  of  the  bubble.  This  is  due  to  the 
change  in  direction  the  supersonic  flow  encounters  when  it  contacts  the  wall. 

Present  Contribution 

The  work  of  Moore  and  Elward  [1992]  was  aimed  at  further  understanding 
the  mechanism  of  shock  formation  near  the  inlet  of  the  tip  clearance  gap.  The 
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flow  structure  was  related  to  the  development  of  compressible  flow  in  sharp-edged 
orifices.  Particular  features  of  interest  include  the  length  scale  of  the  formation 
process  and  the  strength  of  the  shocks  produced.  In  this  technical  note,  the  findings 
of  Moore  and  El  ward  are  summarized. 


Model  of  shock  formation 


Figure  2 shows  the  flow  model  postulated  by  Moore  and  Elward  for  the 
shock  formation.  The  flow  separates  from  the  corner  at  the  tip  gap  inlet. 
Supersonic  freestream  flow  with  a Mach  number  M,  then  overturns  and  intersects 
the  blade  tip  at  a distance  Xj  from  the  corner.  Here  the  flow  is  abruptly  turned 
through  an  angle  5,  giving  rise  to  an  oblique  shock  wave  at  an  angle  a-8  to  the 
surface.  This  shock  formation  process  is  simply  modelled  as  shock  formation  in 
supersonic  flow  over  a wedge.  The  shock  angles  are  then  given  by 
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Effective  wall  location  of  wave  formation 

Figure  3 shows  the  free  streamlines  calculated  for  sharp-edged  orifice  flows 
by  Norwood. 

From  the  free  streamline  results  of  Benson  and  Pool  and  of  Norwood,  Fig. 
4 was  constructed.  This  figure  shows  the  distance,  \t,  from  the  orifice  entrance 
to  the  point  of  intersection  of  the  free  streamline  with  a line  drawn  from  the  orifice 
edge  parallel  to  the  orifice  centerline,  plotted  against  the  freestream  Mach  number. 
This  distance  is  plotted  as  x/w,  or  the  distance  in  orifice  half-widths.  The  figure 
shows  that  at  higher  Mach  numbers,  the  free  streamline  intersects  the  "wall"  closer 
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to  the  orifice  edge.  As  the  Mach  number  decreases,  the  free  streamline  intersects 
the  wall  farther  and  farther  downstream.  The  free  streamline  becomes  parallel  to 
the  wall  at  a pressure  ratio  of  about  0.36  or  M ~ 1.3,  as  suggested  by  the  results 
of  Benson  and  Pool. 

Figure  4 also  plots  the  data  for  five  water  table  cases  against  the  Mach 
number,  M,,  from  Moore  and  Elward.  Two  points  are  plotted  for  each  case,  the 
location  of  the  intersection  of  the  line  of  median  heights  with  the  channel  wall, 
denoted  by  the  symbol  m,  and  the  location  of  the  intersection  of  the  line  of 
maximum  heights  with  the  channel  wall,  denoted  by  the  symbol  p.  The  intersec- 
tions for  the  lines  of  maximum  height  agree  well  with  the  wall  locations  predicted 
from  the  calculated  compressible  flow  free  streamlines.  This  data  lies  in  the  range 

x/w  = 0.9-1. 2.  The  trend  in  the  data  follows  the  predicted  variation  with  free 
stream  Mach  number. 

The  angles  of  turn,  5,  or  equivalently  the  angles  with  which  the  free 
streamlines  intersect  the  wall,  are  plotted  against  freestream  Mach  number  in  Fig. 
5.  The  figure  shows  the  turning  angles  from  Moore  and  Elward  and  the  free 
streamline  angles  of  Benson  and  Pool  and  of  Norwood.  As  was  seen  in  Fig.  4, 
for  Mach  numbers  above  about  1.3,  the  free  streamline  intersects  the  wall.  The 
angle  of  intersection  then  becomes  larger  as  the  Mach  number  increases. 

The  line  of  minimum  Mach  number  required  for  an  attached  shock  for  a 
given  5 is  also  shown  on  Fig.  5.  It  appears  that  the  shock  formation  is  like  that  of 
an  attached  shock  on  a wedge  of  half-angle,  5. 

Another  interpretation  of  the  shock  formation  is  that  it  is  like  turbulent 
reattachment  in  supersonic  flow.  Carriere  [1970]  has  presented  a correlation  of 
experimental  results  for  two-dimensional  flow.  Again  Figure  5 shows  that  this  is 
in  reasonable  agreement  with  the  data,  but  the  trend  is  toward  somewhat  lower 
turning  angles,  or  later  reattachment,  at  higher  Mach  numbers.  This  perhaps 
supports  the  argument  by  Norwood  that  the  phenomenon  of  reattachment  on  his 
inclined  wall  was  a compressibility  effect  rather  than  a frictional  effect. 

Wave  formation  in  compressible  flows 


The  results  in  Figs.  4 and  5 may  be  used  to  predict  oblique  shock  formation 
in  compressible  flows.  For  example,  consider  a flow  with  a free  stream  Mach 
number  of  1.8,  that  is,  a minimum  pressure  p,/p0  = 0.174  (with  k = 1.4). 
Figures  4 and  5 give  the  wave  location  and  turning  angle  as  x/w  = 0.96  and  5 = 
18.8  degrees,  respectively.  Equations  1 and  2 may  then  be  solved  to  get  a = 61.1 

degrees  and  a - 5 = 42.3  degrees.  The  resulting  predicted  flow  is  shown  plotted 
in  Fig.  6. 
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Implications  for  iet  engine  heat  transfer 

The  2-D,  incompressible  turbulent  flow  calculations  of  Moore,  et  al.  [1989] 
showed  an  area  of  enhanced  heat  transfer  on  the  pressure  side  of  a turbine  blade 
tip  as  shown  in  Fig.  7.  The  heat  transfer  was  enhanced  by  up  to  1.8  times  the 
downstream  fully  developed  value  in  the  first  two  to  three  tip  gap  heights.  The 
estimated  intersection  of  the  free  streamline  with  the  sidewall  of  x/w  between  0.4 
and  1.4,  observed  in  Fig.  4,  would  indicate  a shock  forming  within  this  region  of 
already  enhanced  heat  transfer.  The  shock-boundary  layer  interaction  could  serve 
to  further  enhance  the  heat  transfer.  Increased  heat  transfer  would  lead  to  higher 
metal  temperatures  and  increased  rates  of  oxidation  and  material  weight  loss.  This 
would  reduce  both  the  expected  useful  life  of  the  turbine  blade  and  the  turbine 
efficiency. 
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distance  from  orifice  edge  (half  orifice  widths) 


Norwood’s  calculated  free  streamlines  from  a half-slit 
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Figure  4.  Comparison  of  free  streamline  Figure  5.  Compressible  flow  turning  angles 

intersection  with  channel  sidewall 


Figure  6.  Predicted  shock  formation  in  compressible  gas  flow  with  MfrecstreJU11  =1.8 


range  of  locations  of  wave 
formation  in  present  tests 


calculated  with  T — 0 in  blade 

— calculated  with  T«0  at  base  of  blade  tip 
0 *values  measured  by  Metzger  and  Bunker 


Eigure  7,  Enhanced  heat  transfer  to  turbine  blade  tip 
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Axisymmetric  Computational  Fluid  Dynamics  Analysis  of 
Saturn  V/S1-C/F1  Nozzle  and  Plume 

Joseph  H.  Ruf 

NASA/Marshall  Space  Flight  Center 
Huntsville,  AL. 

Abstract 

An  axisymmetric  single  engine  Computational  Fluid 
Dynamics  calculation  of  the  Saturn  V/Sl-C  vehicle  base  region 
and  FI  engine  plume  is  described.  There  were  two  objectives 
of  this  work,  the  first  was  to  calculated  an  axisymmetric 
approximation  of  the  nozzle,  plume  and  base  region  flow  fields 
of  Sl-C/Fl,  relate/scale  this  to  flight  data  and  apply  this  scaling 
factor  to  a NLS/STME  axisymmetric  calculations  from  a parallel 
effort.  The  second  was  to  assess  the  differences  in  FI  and 
STME  plume  shear  layer  development  and  concentration  of 
combustible  gases.  This  second  piece  of  information  was  to  be 
input/supporting  data  for  assumptions  made  in  NLS2  base 
temperature  scaling  methodology  from  which  the  vehicle  base 
thermal  environments  were  being  generated.  The  FI 
calculations  started  at  the  main  combustion  chamber  faceplate 
and  incorporated  the  turbine  exhaust  dump/nozzle  film 
coolant.  The  plume  and  base  region  calculations  were  made  for 
ten  thousand  feet  and  57  thousand  feet  altitude  at  vehicle 
flight  velocity  and  in  stagnant  freestream.  FDNS  was 
implemented  with  a 14  species,  28  reaction  finite  rate 
chemistry  model  plus  a soot  burning  model  for  the  RP-l/LOX 
chemistry.  Nozzle  and  plume  flow  fields  are  shown,  the  plume 
shear  layer  constituents  are  compared  to  a STME  plume. 
Conclusions  are  made  about  the  validity  and  status  of  the 
analysis  and  NLS2  vehicle  base  thermal  environment  definition 
methodology. 
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An  in  house  CFD  effort  to  qualitatively  assess  NLS/STME 
base  heating  rates  was  begun. 

- this  included  similar  axisymmetric  analysis  of 
Saturn  V/S1-C/F1  and  NLS/STME  configurations 
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ase  thermal  environments  were  being  generated. 
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vs.  1 large  and  21  smaller  slots 


FI  MCC  end  Nozzle  brio 
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Chemistry 

- finite  rate,  with  1 4 species  and  28  reactions 

- soot  modeled  as  solid  carbon 


10k ft  Plune  Grid 
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- calculated  thrust  and  Isp 
frozen  flow  +.5% 
finite  rate  +12% 
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Pressure  (psf) 
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ratio  1 .8  to  1 .9  to  1 


tort 
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Appears  to  be  a significant  difference  in  the  plume  shear  layer 
development  between  FI  and  STME  (2nd  objective). 

- Indicates  that  plume  shear  layer  development  and  combustible 
gas  concentrations  are  not  similar,  therefore,  the  NLS2  Base 
Gas  Temperature  Scaling  Methodology  may  be  non  conservative. 
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Computational  Fluid  Dynamic  (CFD)  Analysis  of  Axisymmetric 
Plume  and  Base  Flow  of  a Film/Dump  Cooled  Rocket  Nozzle 


P K.  Tucker 

NASA/Marshall  Space  Flight  Center  (MSFC) 
Marshall  Space  Flight  Center,  AL  35812 

S.  A.  Warsi 

Sverdrup  Technology,  Inc.  (MSFC  Group) 
620  Discovery  Drive 
Huntsville,  AL  35806 


Film/dump  cooling  a rocket  nozzle  with  fuel  rich  gas,  as  in  the  National  Launch  System 
(NLS)  Space  Transportation  Main  Engine  (STME),  adds  potential  complexities  for 
integrating  the  engine  with  the  vehicle.  The  chief  concern  is  that  once  the  film  coolant 
is  exhausted  from  the  nozzle,  conditions  may  exist  during  flight  for  the  fuel-rich  film 
gases  to  be  recirculated  to  the  vehicle  base  region.  The  result  could  be  significantly 
higher  base  temperatures  than  would  be  expected  from  a regeneratively  cooled 
nozzle. 

CFD  analyses  were  conduced  to  augment  classical  scaling  techniques  for  vehicle 
base  environments.  The  FDNS  code  with  finite  rate  chemistry  was  used  to  simulate  a 
single,  axisymmetric  STME  plume  and  the  NLS  base  area.  Parallel  calculations  were 
made  of  the  Saturn  V S-1C/F1  plume  base  area  flows.  The  objective  was  to 
characterize  the  plume/freestream  shear  layer  for  both  vehicles  as  inputs  for  scaling 
the  S-C/F1  flight  data  to  NLS/STME  conditions.  The  code  was  validated  on  high 
speed  flows  with  relevant  physics.  This  paper  contains  the  calculations  for  the 
NLS/STME  plume  for  the  baseline  nozzle  and  a modified  nozzle.  The  modified  nozzle 
was  intended  to  reduce  the  fuel  available  for  recirculation  to  the  vehicle  base  region. 
Plumes  for  both  nozzles  were  calculated  at  lOkFT  and  50kFT. 


CLANK  NOT 


^PAGEji^ 


1457 


INTENTIONALLY  BLANK 


National  Aeronautics  and 
Space  Administration 


UJ 

CO 

< 

m 

08 


o 

DC 

H 

LU 


>- 

CO 

X 

< 


LL 

O 


< ' 
LL 

O 


S CM 
1 8 
c O 


(0 

2 


$ 


1458 


Saif  Wars! 

Sverdrup  Technology,  Inc. 


a> 


co 

O) 

CC 

c 

o 

o> 

0> 

oc 

© 

CO 

CO 

OQ 


5 

UJ 

CC 

UJ 


0> 


o 

J* 

o 


o 

CO 

o 


CL 

< 


o 

<0 

s 

Q. 

E 

o 

O 

C "D 
O CD 
CO 

2.  5 

E ! 

O CD 

O -E 

*->  q 
-2  co 
o co 
IO  CQ 
CO  ^ ;tf 

■=  32  12 

2 O CD 
8 « 
IE  ' ' 


CO 

3 

s 

35 


1459 


CP 


CP 

© 

CC 


O) 

© 

cc 

© 

© 

CO 

CD 

T5 

C 

© 

© 

E 

3 

Q. 

Ill 

2 

h 

CO 

CO 


SB 


« -2 
2 £ 
s £ 

u 

Ji 

« a 
z co 


CO 

LLi 

> 

h 

O 

LLI 

m 

o! 


© 

■o 

3 


© 

O 

c 

o 

o 

c 

3 

H- 

© 

© 

© 

© 

o 

© 

© 

© 

© 

E 


IXJ 

2 

CO 

CO 


© 

© 

© 

c 

© 

C3 


© 

N 

N 

O 

C 

o 

E 


© 

3 

© 

© 

CC 


© 

w 

© 

© 

© 

V. 

+-• 

W 

E 

© 

c 

£ 

<— 

© 

c 

o 

• MM© 

"O 

© 

o 

"O 

© 

CL 

O) 

© 

© 

o 

E 

© 

"O 

3 

c 

© 

"D 

© 

© 

1— 

© 

ID 

© 

c 

© 

*♦— 

M— 

V. 

© 

o 

o 

© 

r~ 

c 

c 

© 

o 

o 

© 

O 

t3 

o 

o 

c 

c 

c 

o 

3 

H- 

3 

o 

© 

'u 

© 

© 

© 

© 

© 

© 

© 

Q. 

< 

< 

CL 

E 

• 

• 

E 

o 

o 

o 

O 

© 

© 


D) 

h— 

c 

1— 

3 

© 

CO 

o> 

■© 

o 

© 

o 

H— 

o 


c 

o 


1460 


Develop  significant  in-house  capability  for 

- Reacting  nozzle  flows 

- Reacting  plumes 

- Complex  base  flows 
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NLS/STME  Plume  and  Base  Region  Results 
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NLS/STME  Plume  and  Base  Region  Results 
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NLS/STME  Plume  and  Base  Region  Results 
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NLS/STME  Plume  and  Base  Region  Results 
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Preliminary  calculations  at  high  altitudes  indicate: 

• at  1/2  radius  downstream  of  nozzle  more  combustibles  (approx.  2:1)  in  NLS/STME 
shear  layer 

• at  1 radius  downstream  of  nozzle  , combustible  ratio  is  about  even 
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NLS  Base  Heating  CFD  Analysis 

Edward,  P.  Ascoli,  Adel  H.  Heiba,  Yann-Fu  Hsu, 
Ronald  R.  Lagnado,  and  Edward  D.  Lynch 

Rockwell  International,  Rocketdyne  Division 


Workshop  for  Computational  Fluid  Dynamic 
Applications  in  Rocket  Propulsion 

April  20-22,  1993 

NASA  Marshall  Space  Flight  Center 


Abstract 

Concerns  raised  over  possible  base  heating  effects  on  the  NLS  1.5  stage  reference 
vehicle  resulted  in  the  use  of  CFD  as  a predictive  analysis  tool.  The  objective 
established  was  to  obtain  good  engineering  solutions  to  describe  the  base  region 
flowfields  at  10,000  ft.  and  50,000  ft.  altitudes.  The  Rockwell  USA  CFD  code  was 
employed  with  a zero-equation  turbulence  model  and  a four  species,  1 step  chemical 
kinetics  package.  Three  solutions  were  generated  for  the  specified  altitudes  on  coarse 
and  fine  grids.  CFD  results  show  the  base  region  flowfields  to  be  highly  three- 
dimensional  in  character.  At  the  10,000  ft.  altitude,  plumes  contract  soon  after  exiting 
the  nozzles  and  do  not  interact  with  each  other.  No  mechanism  was  identified  for 
driving  hot  gas  back  into  the  base  region  and  no  significant  amounts  of  hydrogen  or 
water  were  found  in  the  base  region.  Consequently,  surface  temperatures  were  all 
near  the  ambient  level.  At  50,000  ft.,  the  nozzle  exhaust  plumes  begin  to  interact, 
particularly  those  of  the  two  inboard  engines  which  are  closer  together.  A small 
amount  of  hot  gas  is  recirculated  between  the  inboard  nozzles  near  the  nozzle  exit 
plane.  As  a result,  base  region  surface  temperatures  are  slightly  elevated,  but  still 
remain  well  within  the  design  guideline  of  1000°R. 
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CFD  IDENTIFIED  AS  ALTERNATE  PREDICTIVE  TOOL 
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ALL  CALCULATIONS  DONE  WITH  USA  CODE 


CFD  MODEL  OF  NLS  1.5  STAGE  VEHICLE 
BASE  AND  ENGINE  CLUSTER  (CONT'D) 
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NLS  1.5  STAGE  VEHICLE  BASE  AND  ENGINE  CLUSTER 

MODELED  THROUGH  SYMMETRY 
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1.5  STAGE  REFERENCE  VEHICLE  GRID  TOPOLOGY 
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1.5  STAGE  REFERENCE  VEHICLE  GRID 
TOPOLOGY  (CONT’D) 
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BOUNDARY  CONDITIONS  APPLIED 


I486 


Rockwell  International 

Rocketdyne  Division  cfo  9zo42-o67/o3/slb 


1487 


^ CO 

Z LU 

O g 

t > 
Q O 

z cc 
O CL 
O ' 

d > 

< I— 

5 $ 

o CC 

F w 

< CO 

CD  5 

< o 
5 ° 


< s 

o b 

Q Z 

nr  ul 

O g 

LU  > 

# m 
cc 

< X 

o 3 

o El 

cc  H 

P < 

LL  LU 
LU  I 


LU  ^ 

n “ 

o CO 

2 < 


LL 

o 

LU  CC 

DC  0- 
D Q 
H LU 

a y= 

lu  O 
CL  LU 

S CL 

H DC 

-I  o 


o 

h- 

g 5 
- § 

§ < 
OT  ° 

CC  * 

< 3 

lu  Cl 

I — 
Q < 


w b 

5 LU 


< < LU  CC 

5 2 to  Q 

Q (0 

LU  LL  g 

t h o O 

o o u -1 

LU  CO 
0-  ~ 

C/)  • 


LU  _| 
DC  < 


1488 
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PARTICLE  TRACES  IN  BASE  REGION 

10,000  FOOT  ALTITUDE 
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PARTICLE  TRACES  IN  BASE  REGION  (CONT’D) 

1 0,000  FOOT  ALTITUDE  
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REGION  MACH  CONTOURS  (CONT’D) 

10,000  FOOT  ALTITUDE 
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PLUME  REGION  MACH  CONTOURS 

10,000  FOOT  ALTITUDE 
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PARTICLE  TRACES  FROM  NOZZLE  EDGE 

1 0,000  FOOT  ALTITUDE 


1493 


Rockwell  International 

Rocketdyne  Division 


BASE  REGION  H2  CONTOURS  (CONT’D) 

10,000  FOOT  ALTITUDE 
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REGION  TEMPERATURE  CONTOURS  (CONT'D) 

10,000  FOOT  ALTITUDE 
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BASE  REGION  SURFACE  TEMPERATURE  CONTOURS 

10,000  FOOT  ALTITUDE 

TEMP  (°R) 
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NOZZLE  SURFACES  TEMPERATURES 

10,000  FOOT  ALTITUDE 
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10,000  FOOT  ALTITUDE  RESULT  SUMMARY 

KEY  FLOW  FEATURES  CAPTURED 
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COARSE  AND  FINE  GRID  SOLUTIONS  PREDICT  ESSENTIALLY  SAME 
FLOWS 


BASE  REGION  MACH  CONTOURS  (CONT'D) 

50,000  FOOT  ALTITUDE 
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PARTICLE  TRACES  FROM  NOZZLE  EDGE 

50.000  FOOT  ALTITUDE 


1500 


Rockwell  International 

Rocketdyne  Division  cfd  ga-wz^a^a/sLB 


PARTICLE  TRACES  FROM  NOZZLE  EDGE  (CONT'D) 

50.000  FOOT  ALTITUDE 
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PARTICLE  TRACES  IN  BASE  REGION 

50,000  FOOT  ALTITUDE 
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PARTICLE  TRACES  IN  BASE  REGION  (CONT’D) 

50,000  FOOT  ALTITUDE 
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BASE  REGION  H20  CONTOURS  (CONT'D) 

50,000  FOOT  ALTITUDE 


1504 


Rockwell  International 

Rocketdyne  Division 


BASE  REGION  SURFACE  TEMPERATURE  CONTOURS 

50,000  FOOT  ALTITUDE 

TEMP  (°R) 
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NOZZLE  SURFACES  TEMPERATURES 

50,000  FOOT  ALTITUDE 

OUTER  NOZZLE  CENTER  NOZZLE 
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50,000  FOOT  ALTITUDE  RESULT  SUMMARY 
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CFD  92-042-1 33/D3/SLB 


BASE  FLOW  ENVIRONMENT  DEFINED  FOR  NLS  1.5  STAGE 
REFERENCE  GEOMETRY  AT  10.000  AND  50.000  FT  ALTITUI 
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ABSTRACT: 

CFD  Flowfield  Simulation  of  Delta  Launch  Vehicles  in  a Power-on  Configuration  o 

r 

Workshop  for  Computational  Fluid  Dynamic  (CFD)  Applications  in  Rocket  Propulsion, 
George  C.  Marshall  Space  Flight  Center,  April  20-22, 1993 

By:  D.  L.  Pavish+,  T.  P.  Gielda*,  B.  K.  Soni  **,  J.  E.  Deese++,  and  R.  K.  Agarwal++ 


s^-3i 


This  paper  summarizes  recent  work  at  McDonnell  Douglas  Aerospace  (MDA)  to  develop 
and  validate  computational  fluid  dynamic  (CFD)  simulations  of  under  expanded  rocket 
plume  external  flowfields  for  multibody  expendable  launch  vehicles  (ELVs).  Multi- 
engine  reacting  gas  flowfield  predictions  of  ELV  base  pressures  are  needed  to  define 
vehicle  base  drag  and  base  heating  rates  for  sizing  external  nozzle  and  base  region 
insulation  thicknesses.  Previous  ELV  design  programs  used  expensive  multibody 
power-on  wind  tunnel  tests  that  employed  chamber/nozzle  injected  high  pressure  cold 
or  hot-air.  Base  heating  and  pressure  measurements  were  belatedly  made  during  the 
first  flights  of  past  ELVs  to  correct  estimates  from  semi-empirical  engineering  models  or 
scale  model  tests. 


Presently,  CFD  methods  for  use  in  ELV  design  are  being  jointly  developed  at  the  Space 
Transportation  Division  (MDA-STD)  and  New  Aircraft  Missiles  Divison  (MDA-NAMD). 
An  explicit  three  dimensional,  zonal,  finite-volume,  full  Navier-Stokes  (FNS)  solver  with 
finite  rate  hydrocarbon/air  and  aluminum  combustion  kinetics  was  developed  to 
accurately  compute  ELV  power-on  flowfields.  Mississippi  State  University’s  GENIE++ 
general  purpose  interactive  grid  generation  code  was  chosen  to  create  zonal,  finite 
volume  viscous  grids.  Axisymmetric,  time  dependent,  turbulent  CFD  simulations  of  a 
Delta  DSV-2A  vehicle  with  a MB-3  liquid  main  engine  burning  RJ-1/LOX  were  first 

completed.  Hydrocabon  chemical  kinetics  and  a k-e  turbulence  model  were  employed 
and  predictions  were  validated  with  flight  measurements  of  base  pressure  and 
temperature.  Zonal  internal/external  grids  were  created  for  a Delta  DSV-2C  vehicle  with 
a MB-3  and  three  Castor-1  solid  motors  burning  and  a Delta-ll  with  an  RS-27  main 
engine  (LOX/RP-1)  and  9 GEMs  attached/6  burning.  Cold  air,  time  dependent  FNS 
calculations  were  performed  for  DSV-2C  during  1992.  Single  phase  simulations  that 
employ  finite  rate  hyrocarbon  and  aluminum  (solid  fuel)  combustion  chemistry  are 
currently  in  progress.  Reliable  and  efficient  Eulerian  algorithms  are  needed  to  model 
two  phase  (solid-gas)  momentum  and  energy  transfer  mechanisms  for  solid  motor  fuel 
combustion  products. 


+ McDonnell  Douglas  Aerospace-Space  Transportation  Division,  Huntingon  Beach, Ca. 

++  McDonnell  Douglas  Aerospace-New  Aircraft  Missile  Products  Division,  St.  Louis,  Mo. 
* Presently  at  Ford  Motor  Co.,  Detroit,  Mi. 

**  Mississippi  State  Univ.  NSF/Eng.  Research  Center,  Strarkville,  Miss. 
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Finite  rate  combustion  chemistry  routines  (explicit) 

hydrocarbon  (RJ-1  or  RP-1)  and  aluminum  based  solid  motor  fuel 

Robust,  well  validated  CFD  code  for  ELV  power-off  flowfields 
( see  AIAA-91-1727  and  AIAA-92-2681 ) 


APPROACH  (Continued) 


Assume  single  phase  gas  reactions  w/o  surface  combustion 

Utilize  Delta  launch  vehicle  flight  measurements  of  base 
surface  pressures,  temperatures,  heating  rates 
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DELTA  DSV  FLIGHT  419  RECONSTRUCTED  TRAJECTORY 
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Velocity  (kft/sec) 
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Altitude  = 160  kft, 
Mach  Number  = 5.37 
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Altitude  = 160  kft, 
Mach  Number  = 5.37 


DSV  FLOWFIELD  C02  MASS  FRACTIONS 

MDA-Space  Transportation  Division 
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Altitude  = 160  kft, 
Mach  Number  = 5.37 


DELTA  DSV  BASE  SURFACE  TEMPERATURE 
COMPARISON  OF  CFD  PREDICTION  VS  FLIGHT  MEASUREMENTS 
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Mach  Number  = 1.2  Mach  Number  =2.2 
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Employ  aluminum  solid  fuel  finite  rate  combustion  chemistry 


SUMMARY  AND  CONCLUSIONS 
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NUMERICAL  STUDY  OP  THE  SSME  NOZZLE  FLOW  PIELDS 
DURING  TRANSIENT  OPERATIONS  - A COMPARISON  OF  THE  ANIMATED  y,  ^ 

RESULTS  WITH  TEST 


Ten-See  Wang 

Computational  Fluid  Dynamics  Branch 
NASA  - Marshall  Space  Flight  Center 
Marshall  Space  Flight  Center,  AL  35812 


Catherine  Dumas 
SVERDRUP  Technology,  Inc. 
Huntsville,  AL  35806 


Abstract 


A computational  fluid  dynamics  (CFD)  model  has  been  applied 
to  study  the  transient  flow  phenomena  of  the  nozzle  and  exhaust 
plume  of  the  Space  Shuttle  Main  Engine  (SSME),  fired  at  sea  level. 
The  CFD  model  is  a time  accurate,  pressure  based,  reactive  flow 
solver.  A six— species  hydrogen/oxygen  equilibrium  chemistry  is 
used  to  describe  the  chemical-thermodynamics.  An  adaptive 
upwinding  scheme  is  employed  for  the  spatial  discretization,  and 
® P^s^ictor , multiple  corrector  method  is  used  for  the  temporal 
solution.  Both  engine  start— up  and  shut-down  processes  were 
simulated.  The  elapse  time  is  approximately  five  seconds  for  both 
cases.  The  computed  results  were  animated  and  compared  with  the 
test.  The  images  for  the  animation  were  created  with  PLOT3D  and 
FAST  and  then  animated  with  ABEKAS.  The  hysteresis  effects,  and 
the  issues  of  free— shock  separation,  restricted— shock  separation 
and  the  end— effects  were  addressed. 
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Numerical  Study  of  the  SSME  Nozzle  Flowfields 

during  Transient  Operations 
- A Comparison  of  Animated  Results  with  Test 
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Open  Forum:  Combustion  - Nozzles/plumes  - Comparison 

MSFC,  Alabama 
April  22, 1993 
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• PL0T3D  and  FAST  generated  images 

• ABEKAS  generated  animatioins 
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• A comparison  of  computed  thrust  coefficients  with  those  of  measurement 

- Start-up  transient 

- Shut-down  transient 
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N95- 23435 
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Enhanced  design  / analysis  tools  available  for  future  application 


Volute  Design  Process  Outline 


Advanced  Concept  Volute  Geometry 


Exhaust  Volute  Manifold  Design  Requirements 
High  Turbine  Stage  Loading  Creates  Challenging  Goals 
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Turbin*  loading  coefficient 


Meanline  Analysis 

First  Pass  Estimate  of  Geometry  and  Performance 
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Parametric  Studies 

- Performance  Optimization  and  Sensitivity 


Geometry  and  Mesh  Generation 

Rules  Based  Design  Program  Used  To  Create  Volute  Geometry 
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Baseline  Volute  3D  Flow  Analysis 
Circumferential  Static  Pressure  Distribution 
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Circ.  Position  (deg.) 
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Circ.  Position  (deg.) 


Baseline  Volute 
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Integrated  CFD  Modeling  of  Components 

Computational  Mesh  for  Inlet  Volute,  Turbine  Stage  and  Exhaust  Manifold 
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Time  Accurate  3D  Turbine  Stage  Flow  Analysis 
Euler  (w  / shear),  20  Vanes,  and  42  Blades 
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Integrated  Time  Accurate  CFD  Modeling  of  Components 

Turbine  Stage  (20  vanes,  42  Blades)  and  Exhaust  Manifold  (42  Elements) 
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Circ.  Position  (deg.) 


Time  Accurate  3D  Turbine  Stage  Flow  Analvsi 
Euler  (w  / shear),  20  Vanes,  and  42  Blades 
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Circ.  Position  (deg.) 


Integrated  Time  Accurate  CFD  Modeling  of  Components 
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Circ.  Position  (deg.) 


Integrated  Time  Accurate  CFD  Modeling  of  Components 

Turbine  Stage  (20  Vanes,  42  Blades)  and  Exhaust  Manifold  (42  Elements) 
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Circ.  Position  (deg.) 


Integrated  Time  Accurate  CFD  Modeling  of  Components 
Turbine  Stage  (20  Vanes,  42  Blades)  and  Exhaust  Manifold  (42  Elements) 
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Aero  Design  and  Analysis  of  a Highly  Loaded  Turbine  Exhaust  Volute  Manifold 

SUMMARY 
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Results  To  Be  Utilized  In  Advanced  Concept  Definition 
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CFD  Analysis  of  Turbopump  Volutes 

Edward,  P.  Ascoli,  Daniel  C.  Chan,  Armen  Darian, 
Wayne  W.  Hsu,  and  Ken  Tran 

Rockwell  International,  Rocketdyne  Division 


Workshop  for  Computational  Fluid  Dynamic 
Applications  in  Rocket  Propulsion 

April  20-22,  1993 

NASA  Marshall  Space  Flight  Center 
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Abstract 

An  effort  is  underway  to  develop  a procedure  for  the  regular  use  of  CFD  analysis  in  the 
design  of  turbopump  volutes.  Airflow  data  to  be  taken  at  NASA  Marshall  will  be  used 
to  validate  the  CFD  code  and  overall  procedure.  Initial  focus  has  been  on 
preprocessing  (geometry  creation,  translation,  and  grid  generation).  Volute 
geometries  have  been  acquired  electronically  and  imported  into  the  CATIA  CAD 
system  and  RAGGS  (Rockwell  Automated  Grid  Generation  System)  via  the  IGES 
standard.  An  initial  grid  topology  has  been  identified  and  grids  have  been  constructed 
for  turbine  inlet  and  discharge  volutes.  For  CFD  analysis  of  volutes  to  be  used 
regularly,  a procedure  must  be  defined  to  meet  engineering  design  needs  in  a timely 
manner.  Thus,  a compromise  must  be  established  between  making  geometric 
approximations,  the  selection  of  grid  topologies,  and  possible  CFD  code 
enhancements.  While  the  initial  grid  developed  approximated  the  volute  tongue  with  a 
zero  thickness,  final  computations  should  more  accurately  account  for  the  geometry  in 
this  region.  Additionally,  grid  topologies  will  be  explored  to  minimize  skewness  and 
high  aspect  ratio  cells  that  can  affect  solution  accuracy  and  slow  code  convergence. 
Finally,  as  appropriate,  code  modifications  will  be  made  to  allow  for  new  grid 
topologies  in  an  effort  to  expedite  the  overall  CFD  analysis  process. 
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CFD  ANALYSIS  OF  TURBOPUMP  VOLUTES 
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April  20-22, 1993 

NASA  Marshall  Space  Flight  Center 


TASK  OBJECTIVES 
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PREPROCESSING  PROCEDURE  ESTABLISHED 


1559 


CFD  93-013(C)-004/Dl/SLB 


ELECTRONIC  DESIGN  & GRID  GENERATION  TOOLS 
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INTERACTIVE  USER  INTERFACE 
ACCEPTS  IGES  GEOMETRY  FILES 


EXTENSIVE  TRANSLATION/VERIFICATION  REQUIRED 
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THREE  ZONE  GRID  DEVELOPED 


1566 


Rockwell  International 

Rocketdyne  Division  cfd93-oi3<c)-oo8/oi/slb 


THREE-ZONE  VOLUTE  TOPOLOGY 
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ZONE  1 - INLET  GRID 
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Rocketdyne  Division  CFD  92-059-01 6/D4/MMS 


ZONE  2 - TAPERED  MANIFOLD 
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Z0NE.3  - ANNUAL  DISCHARGE 
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SAME  TOPOLOGY  USED  FOR  EXIT  VOLUTE 
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EXIT  VOLUTE  ZONE  1 
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DESIGN  PROCEDURE  MUST  MEET  ENGINEERING 
ACCURACY,  CYCLE  TIME  REQUIREMENTS 
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OVERALL  PROCEDURE 


GEOMETRIC  APPROXIMATIONS 
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GRID  "SMEARING"  POSSIBLE  WITH  SOME  TOPOLOGIES 
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SUMMARY 
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Three-Dimensional  Viscous  Flow  Analysis 
Inside  a Turbine  Volute 


C.Hah,  J.  Loellbach,  and  D.  A.  Green wald 
NASA  Lewis  Research  Center 

L.  Griffin  and  J.  Ruf 
NASA  Marshall  Space  Flight  Center 


A three-dimensional  numerical  method  has  been  developed 
to  analyze  the  complex  flow  field  inside  a turbine  volute. 

Comparisons  are  made  between  solutions  with  different  boundary 
conditions. 
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Consortium  Turbine  and  Exit  Volute 
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Approach 
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TIP  GRID 
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MACH-NUMBER  AND  VELOCITY  VECTORS  (MID-SPAN) 
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Surface  Pressure  Contour  on  Oxidizer  Turbine  Rotor 


Velocity  Vectors  on  Oxidizer  Turbine  Rotor 


Consortium  Turbine  and  Exit  Volute 
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Computational  grid 
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Surface  Triangulation  of  Oxidizer  Turbine  Volute 


BOUNDARY  CONDITIONS 
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HAH3D  Viscous  Flow  Code 
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Velocity  Vectors  Near  Inner  Wall 
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Velocity  Vectors  Near  Outer  Wall 


Crossflow  Velocity  Vectors 
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Surface  Static  Pressure 
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Velocity  Vectors  on  Oxidizer  Turbine  Volute 


1604 


Velocity  Vectors  on  Oxidizer  Turbine  Volute 
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Surface  Pressure  ContotifWrfOxidizer  Turbine  Volute 


Observation  and  Future  Study 


’O 

2 

m 

* 

© 

E 

u 

a> 

3 

w 

C0 

I 

£ 

Q 


a 


CO 

CO 

0> 

o 

CJ 


2 

!>2 

O 


0> 

u 

3 

w 

■o 

a 

2 

CO 

T3 

a> 

u 

£ 

2 

o 

W 

E 

2 

u 

s 

■*-> 

m 

>% 

U 

a> 

u 

a> 

a «« 
a>  ^ 

CO 

C/3 

C/3 

§S 

a>  o 

£ -c 

a> 

u 

S 

■w 

C-> 

« 2 
■°  s 

a> 

z 

c 

c 

wo 

■5 

<D 

U 

*&  2 
5 r’S 

*3 

2 

•33 

CU 

C-> 

C2 

Q 

u 

jy 

a>  "2 
!>  « 

”u 

a> 

o> 

J3 

© 

> 

*3 

w 

<3 

L« 

-*->  -5 

2 2 

2 3 

a C 

J> 

u 

a> 

JS 

w 

u 

cn 

i 

2 co 
O 2 

u 

2 

e£ 

2 

U P 

o 

o 

O 

1606 


N95- 23438 


WORKSHOP  FOR  CFD  APPLICATIONS  IN  ROCKET  PROPULSION 
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PHASE  II HGM  AIR  FLOW  TESTS  IN  SUPPORT  OF  HEX  VANE  INVESTIGATION 

G.  B.  Cox,  Jr. 

L.  L.  Steele 
D.  W.  Eisenhart 

Pratt  & Whitney/Government  Engines  & Space  Propulsion 
West  Palm  Beach,  Fla. 

Following  the  start  of  SSME  certification  testing  for  the  Pratt  & Whitney  Alternate  Turbopump  Development  (ATD) 
High  Pressure  Oxidizer  Turbopump  (HPOTP),  cracking  of  the  leading  edge  of  the  inner  HEX  vane  was 
experienced.  The  HEX  vane,  at  the  inlet  of  the  oxidizer  bowi  in  the  Hot  Gas  Manifold  (HGM),  accepts  the  HPOTP 
turbine  discharge  flow  and  turns  it  toward  the  Gaseous  Oxidizer  Heat  Exchanger  (GOX  HEX)  coil.  The  cracking 
consistently  initiated  over  a specific  circumferential  region  of  the  hex  vane,  with  other  circumferential  locations 
appearing  with  increased  run  time.  Since  cracking  had  not  to  date  been  seen  with  the  baseline  HPOTP,  a fluid- 
structural  interaction  involving  the  ATD  HPOTP  turbine  exit  flowfield  and  the  HEX  inner  vane  was  suspected. 

As  part  of  NASA  contract  NAS8-36801 , Pratt  & Whitney  conducted  air  flow  tests  of  the  ATD  HPOTP  turbine 
turnaround  duct  flowpath  in  the  MSFC  Phase  II  HGM  air  flow  model.  These  tests  included  HEX  vane  strain  gages 
and  additional  fluctuating  pressure  gages  in  the  turnaround  duct  and  HEX  vane  flowpath  area.  Three-dimensional 
flow  probe  measurements  at  two  stations  downstream  of  the  turbine  simulator  exit  plane  were  also  made. 
Modifications  to  the  HPOTP  turbine  simulator  investigated  the  effects  on  turbine  exit  flow  profile  and  velocity 
components,  with  the  objective  of  reproducing  flow  conditions  calculated  for  the  actual  ATD  HPOTP  hardware. 
Testing  was  done  at  the  MSFC  SSME  Dynamic  Fluid  Air  Flow  (Dual-Leg)  Facility,  at  air  supply  pressures  between 
50  and  250  psia.  Combinations  of  turbine  exit  Mach  number  and  pressure  level  were  run  to  investigate  the  effect 
of  flow  regime. 

Information  presented  includes 

1)  Descriptions  of  turbine  simulator  modifications  to  produce  the  desired  flow  environment. 

2)  T ypes  and  locations  for  instrumentation  added  to  the  flow  model  for  improved  diagnostic 
capability. 

3)  Evaluation  of  the  effect  of  changes  to  the  turbine  simulator  flowpath  on  the  turbine  exit  flow 
environment. 

4)  Comparison  of  the  experimental  turbine  exit  flow  environment  to  the  environment  calculated  for 
the  ATD  HPOTP. 
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SSME  ALTERNATE  TURBOPUMP 
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3-DIMENSIONAL  FLOW  PROBES 
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MATERIAL  BETWEEN  SWIRL 
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INITIAL  ORIFICE  PLATE  REVISIONS  RADIAL  PROFILES 

Significant  Variation  In  Mass  Flow  Profiles  Obtained 
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REVISED  HPOTP  TURBINE  SIMULATOR  FLOW  PROFILES 
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NON-INTRUSIVE  MEASUREMENTS  IN  A ROCKET 
ENGINE  COMBUSTOR 


S.  Farhangi,  V.  T.  Gylys,  R.  J.  Jensen 
Rockwell  International/Rocketdyne  Div. 
Canoga  Park,  California  91303 


ABSTRACT 

In  recent  years  analytical  tools  to  characterize  combustor  flow  have 
been  developed  in  order  to  support  design.  To  facilitate  anchoring  of 
combustion  related  physical  models  and  the  CFD  codes  in  which  they  are 
incorporated  considerable  development  and  application  of  non-intrusive 
combustion  diagnostic  capabilities  has  occurred.  Raman  spectroscopy  can 
be  used  to  simultaneously  detect  all  polyatomic  molecules  present  in 
significant  concentrations  and  to  determine  gas  temperature.  This  is 
because  all  molecules  possess  a distinct  temperature  dependent  Raman 
spectrum . 

A multi-point  diagnostic  system  for  non-intrusive  temperature  and 
species  profiling  in  rocket  engines  has  been  developed  at  Rocketdyne. 
In  the  present  effort,  the  system  has  been  undergoing  validation  for 
application  to  rocket  engine  component  testing.  A 4 inch  diameter 
windowed  combustor  with  a coaxial  gas-gas  injector  was  chosen  for  this 
series  of  validation  experiments.  Initially  an  excimer-pumped  tunable 
dye  laser  and  later  a solid  state  Nd-Yag  laser  served  as  excitation 
sources.  The  Raman  signal  was  dispersed  by  a monochromator  and 
detected  by  a gated,  intensified  Charged  Coupled  Device  (CCD)  array. 

Experiments  were  carried  out  prior  to  each  series  of  hot  fire  tests  to 
ensure  that  the  Raman  signal  detected  was  due  to  a spontaneous  rather 
than  a stimulated  Raman  emission  process.  Over  sixty  hot  fire  tests 
were  conducted  during  the  first  series  of  tests  with  the  excimer/dye 
laser.  All  hot  fire  testing  was  at  a mixture  ratio  of  0.5  and  chamber 
pressures  of  -100  and  -300  psia.  The  Raman  spectra  of  hydrogen,  water 
vapor  and  oxygen  recorded  during  single  element  hot  fire  tests  were 
reduced  and  analyzed.  A significant  achievement  was  the  attainment  of 
single  shot  Raman  spectra  in  cold  flow  tests.  Unfortunately,  the  single 
shot  signal-to-noise  ratio  detoriated  to  an  unacceptable  level  during 
the  hot  fire  testing.  Attempts  to  obtain  temperature  data  from  the 
hydrogen  Ql-branch  profiles  obtained  in  hot  fire  tests  suggest  that 
potentially  complicating  factors  may  render  the  approach  of  averaging 
data  on  the  photodiode  array  invalid.  A second  series  of  hot  fire  tests 
was  conducted  with  a 4 element  coaxial  injector  using  the  Nd-Yag  laser. 
A very  compact  and  portable  diagnostics  set  up  was  assembled  for  ease  of 
alignment,  relocation  and  flexibility.  Measurements  were  made  at 
several  regions  in  the  chamber  in  order  to  map  concentration  profiles. 
High  spatial  resolution  and  improved  signal  to  noise  characteristics 
were  demonstrated. 
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VARIATION  OF  SPECIES  MOLE  FRACTION  WITH  RADIAL  DISRANCE 
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VARIATION  OF  SPECIES  MOLE  FRACTION  WITH  RADIAL  DISRANCE 
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Use  of  Computational  Fluid  Dynamics  (CFD)  codes,  prevalent  in  the  rocket  engine 
turbomachinery  industry,  necessitates  data  of  sufficient  quality  and  quantity  to  benchmark  computational 
codes.  Existing  data  bases  for  typical  rocket  engine  configurations,  in  particular  impellers,  are  limited. 
In  addition,  traditional  data  acquisition  methods  have  several  limitations:  typically  transducer 
uncertainties  are  0.5  percent  of  transducer  full  scale  and  traditional  pressure  probes  are  unable  to 
provide  flow  characteristics  in  the  circumferential  (blade-to-blade)  direction.  Laser  velocimetry 
circumvents  these  limitations  by  providing  + 0.5  percent  uncertainty  in  flow  velocity  and  + 0.5  degree 
uncertainty  in  flow  angle.  The  percent  of  uncertainty  in  flow  velocity  is  based  on  the  measured  value,  not 
full  range  capability.  The  laser  electronics  multiple  partitioning  capability  allows  data  acquired  between 
blades  as  the  impeller  rotates,  to  be  analyzed  separately,  thus  providing  blade-to-blade  flow 
characterization.  Unlike  some  probes,  the  non-intrusive  measurements  made  with  the  laser  velocimeter 
does  not  disturb  the  flow. 

To  this  end1(  and  under  Contract  (NAS8-38864)  to  the  National  Aeronautics  and  Space 
Administration  (NASA)  at  Marshall  Space  Flight  Center  (MSFC),  an  extensive  test  program  was 
undertaken  at  Rocketdyne.  Impellers  from  two  different  generic  rocket  engine  pump  configurations  were 
examined.  The  impellers  represent  different  spectrums  of  pump  design:  the  Space  Shuttle  Main  Engine 
(SSME)  high  pressure  fuel  turbopump  (HPFTP)  impeller  was  designed  in  the  1970*5,  the  Consortium  for 
CFD  Application  in  Propulsion  Technology  Pump  Stage  Technology  Team  (Pump  Consortium)  optimized 
impeller  was  designed  with  the  aid  of  modern  computing  techniques.  The  tester  configuration  for  each  of 
the  impellers  consisted  of  an  axial  inlet,  an  inducer,  a diffuser,  and  a crossover  discharge. 

While  the  tested  configurations  were  carefully  chosen  to  be  representative  of  generic  rocket 
engine  pumps,  several  features  of  both  testers  were  intentionally  atypical.  A crossover  discharge, 
downstream  of  the  impeller,  rather  than  a volute  discharge  was  used  to  minimize  asymmetric  flow 
conditions  that  might  be  reflected  in  the  impeller  discharge  flow  data.  Impeller  shroud  wear  ring  radial 
clearances  were  purposely  close  to  minimize  leakage  flow,  thus  increasing  confidence  in  using  the  inlet 
data  as  an  input  to  CFD  programs. 

The  empirical  study  extensively  examined  the  flow  fields  of  the  two  impellers  via  performance  of 
laser  two-focus  velocimeter  surveys  in  an  axial  plane  upstream  of  the  impellers  and  in  multiple  radial 
planes  downstream  of  the  impellers.  Both  studies  were  performed  at  the  impeller  design  flow 
coefficients. 


Inlet  laser  surveys  that  provide  CFD  code  inlet  boundary  conditions  were  performed  in  one  axial 
plane,  with  ten  radial  locations  surveyed.  Three  wall  static  pressures,  positioned  circumferentially  around 
the  impeller  inlet  were  used  to  identify  asymmetrical  pressure  distributions  in  the  inlet  survey  plane.  The 
impeller  discharge  flow  characterization  consisted  of  three  radial  planes  for  the  SSME  HPFTP  impeller 
and  two  radial  planes  for  the  Pump  Consortium  optimized  impeller.  Housing  wall  static  pressures  were 
placed  to  correspond  to  the  radial  locations  surveyed  with  the  laser  velocimeter.  Between  five  and 
thirteen  axial  stations  across  the  discharge  channel  width  were  examined  in  each  radial  plane  during  the 
extensive  flow  mapping. 

The  largely  successful  empirical  flow  characterization  of  two  different  impellers  resulted  in  a 
substantial  contribution  to  the  limited  existing  data  base,  and  yielded  accurate  data  for  CFD  code 
benchmarking. 
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PUMP  CFD  CODE  VALIDATION  TESTS 
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Two  Pump  Configurations  Tested 
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PUMP  CFD  CODE  VALIDATION  TESTS 
Configurations  1 & 2 Laser  Two-Focus  Velocimeter  Surveys 
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LASER  VELOCIMETER  DATA  SURVEY  LOCATIONS 


PUMP  CFD  CODE  VALIDATION  TESTS 
Consortium  Baseline  Impeller 
Laser  Velocimeter  Data  Survey  Locations 
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PUMP  CFD  CODE  VALIDATION  TESTS 
Configurations  1 & 2 Laser  Two-Focus  Velocimeter  Surveys 
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Configuration  2 - two  radial  survey  planes 

• immediately  downstream  of  impeller  discharge  - 0.51 38  nondimensional  radius 

• Diffuser-side  discharge  plane  - 0.5507  nondimensional  radius 
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SSME  HPFTP  IMPELLER  DISCHARGE 
LASER  VELOCIMETER  DATA  CIRCUMFERENTIAL  REFERENCE 
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CONSORTIUM  BASELINE  IMPELLER  DISCHARGE 
LASER  VELOCIMETER  DATA  CIRCUMFERENTIAL  POSITIONS 
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PUMP  CFD  CODE  VALIDATION  TESTS 
Configuration  1 & 2 Laser  Two-Focus  Velocimeter  Surveys 
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PUMP  CFD  CODE  VALIDATION  TESTS 
Configuration  1 - Laser  Velocimeter  Survey 
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PUMP  CFD  CODE  VALIDATION  TESTS 
SSME  HPFTP  Impeller  Inlet  Laser  Survey 

Impeller  Inlet  Laser  Survey  Plane 
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Impeller  Inlet  Flow  Coefficient 


PUMP  CFD  CODE  VALIDATION  TESTS 
SSME  HPFTP  Impeller  Inlet  Laser  Survey 

Impeller  Discharge  Laser  Survey  Plane  2 (Middle) 
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PUMP  CFD  CODE  VALIDATION  TESTS 
SSME  HPFTP  Impeller  Discharge  Laser  Survey  Plane  1 (Impeller) 
Test  Impeller  Inlet  Flow  Coefficient  = 0.256  

eller  Inlet  Design  Flow  Coefficient  = 0.256  t| 


Data  Windows  in  60  Degree  Arc 


PUMP  CFD  CODE  VALIDATION  TESTS 
SSME  HPFTP  Impeller  Discharge  Laser  Survey  Plane  1 (Impeller) 
Test  Impeller  Inlet  Flow  Coefficient  = 0.256 
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Data  Windows  in  60  Degree  Arc 
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Percent  Impeller  Inlet  Annulus  Height 
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Rockwell  International  Percent  Impeller  Inlet  Annulus  Height 
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PUMP  CFD  CODE  VALIDATION  TESTS 
Configuration  2 - Laser  Velocimeter  Survey 
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PUMP  CFD  CODE  VALIDATION  TESTS 

Consortium  Baseline  Imp.  Inlet  Survey 

Impeller  Inlet  Laser  Survey  Plane 
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0.10  0.12  0.14  0.16 

Impeller  Inlet  Flow  Coefficient 


PUMP  CFD  CODE  VALIDATION  TESTS 
Consortium  Baseline  Imp.  Inlet  Survey 
Impeller  Discharge  Laser  Survey  Plane  2 (Middle) 
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PUMP  CFD  CODE  VALIDATION  TESTS 

Consortium  Baseline  Imp.  Inlet  Survey 

Test  Impeller  Inlet  Flow  Coefficient  = 0.144 
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PUMP  CFD  CODE  VALIDATION  TESTS 
Consortium  Baseline  Imp.  Inlet  Survey 
Test  Impeller  Inlet  Flow  Coefficient  = 0.144 
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PUMP  CFD  CODE  VALIDATION  TESTS 
Consortium  Baseline  Imp.  Discharge  Survey  (Impeller) 
Test  Impeller  Inlet  Flow  Coefficient  = 0.144 
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PUMP  CFD  CODE  VALIDATION  TESTS 

Consortium  Baseline  Imp.  Inlet  Survey  Nondlmenslonal  Axial 
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PUMP  CFD  CODE  VALIDATION  TESTS 
Consortium  Baseline  Impeller  Discharge  Laser  Survey 

Nondimensional  Radial  Planes  0.5138, 0.5597,  Impeller  Inlet  Flow  Coefficient  = 0.144 
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with  Smooth  Rotor  Blades 


N95-  23441 


Susan  T.  Hudson 

NASA-Mar shall  Space  Flight  Center 
ED3 4 /Experimental  Fluid  Dynamics  Branch 
Huntsville,  AL  35812 


s&- 37 

ms  1 1 oca  i 


Abstract 

Several  tests  of  the  Rocketdyne  configuration  of  the  Space 
Shuttle  Main  Engine  (SSME)  High  Pressure  Fuel  Turbopump  (HPFTP) 
Turbine  have  been  completed  in  the  Turbine  Test  Equipment  (TTE)  at 
Marshall  Space  Flight  Center.  The  tests  involved  using  scaled 
performance  parameters  and  model  measurements  to  predict  the 
performance  of  the  turbine.  The  overall  performance  has  been  the 
primary  objective  of  the  tests  to  date,  but  more  detailed 
measurements  are  also  of  interest.  During  the  most  recent  test  of 
the  Rocketdyne  configuration  of  the  HPFTP  turbine  with  smooth  rotor 
blades,  several  different  measurement  techniques  were  used  to  study 
the  turbine  inlet  and  exit  velocity  profiles,  boundary  layer 
thicknesses,  turbulence  intensities,  etc.  Data  has  been  obtained 
using  various  hot  film  probes  and  three-hole  cobra  probes.  Laser 
Velocimeter  measurements  were  also  made.  The  test  plan  and  test 
data  will  be  presented  and  discussed  as  well  as  lessons  learned  on 
how  to  obtain  the  various  types  of  data. 
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Outline 
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• Conclusions 


Test  Objectives 
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• Determine  effect  of  Reynolds  number  on  measurements 


TTE  Facts 
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Pressure  Systems,  Inc.  (PSI)  System 
Hewlett-Packard  (HP)  3852S  data  acquisition  system 
VAX  — database  and  data  analysis 
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FulfScale  SSME  FuelTurbine  TesTWodel 


Phase  IA  Confisurati 
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COMPARISON  OF  COATED,  UNCOATED,  & POLISHED 

SSME  HPFTP  TURBINE  BLADES 
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Phase  IA  Instrumentation 
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Spatial  Resolution  of  Measurements 
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LASER  VELOCIMETER  SYSTEM  SCHEMATIC 
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SEEGMILLER  ATOMIZER 
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SSME  HPFT  INLET  VELOCITY  PROFILE  AND  TURBULENCE  INTENSITY 

(104%  POWER  LEVEL-HOT  FILM-CYLINDER) 
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SSME  HPFT  EXIT  VELOCITY  PROFILE  AND  TURBULENCE  INTENSITY 

(104°/.  POWER  LEVEL-HOT  FILM-CYLINDER) 
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SSME  HPFT  EXIT  VELOCITY  PROFILES  (104%  POWER  LE 
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SSME  HPFT  Exit  Swirl  Versus  Span  and  Circumferential  Position 

(104%  Power  Level) 


(saajBap)  |j;ms 
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Note:  Positive  swirl  is  clockwise  looking  from  turbine  inlet  to  exit. 


SSME  HPFT  Exit  Velocity  Versus  Span  and  Circumferential  Position 

(104%  Power  Level) 


(sdj)  AjpoiaA  amiosqv 
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Development  of  an  Algebraic  Stress/Two-Layer  Model 
for  Calculating  Thrust  Chamber  Flow  Fields 


By 

C.  P.  Chen,  H.  M.  Shang  and  J.  Huang 

Department  of  Chemical  Engineering 
University  of  Alabama  in  Huntsville 
Huntsville,  AL  35899 


ABSTRACT 
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Following  the  consensus  of  a workshop  in  Turbulence  Modeling  for  Liquid  Rocket 
Thrust  Chambers,  the  current  effort  was  undertaken  to  study  the  effects  of  second-order 
closure  on  the  predictions  of  thermochemical  flow  fields.  To  reduce  the  instability  and 
computational  intensity  of  the  full  second-order  Reynolds  Stress  Model,  an  Algebraic 
Stress  Model  (ASM)  coupled  with  a two-layer  near  wall  treatment  was  developed.  Various 
test  problems,  including  the  compressible  boundary  layer  with  adiabatic  and  cooled  walls, 
recirculating  flows,  swirling  flows  and  the  entire  SSME  nozzle  flow  were  studied  to  assess 
the  performance  of  the  current  model.  Detailed  calculations  for  the  SSME  exit  wall  flow 
around  the  nozzle  manifold  were  executed.  As  to  the  overall  flow  predictions,  the  ASM 
removes  another  assumption  for  appropriate  comparison  with  experimental  data,  to  account 
for  the  non-isotropic  turbulence  effects. 
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APPROACH  — Pressure-Strain  Term 
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two-layer  model 
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Two-Layer  Wall  Treatment 
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METHODOLOGIES 
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Time  Marching 


Validations 
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Compressible  Flat  Plate  Flow 
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Fig.  1 Variation  of  Cf/Cf0  with  Moo  for  adiabatic  wall  boundary  condition. 
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2 Variati 


Compressible  Flat  Plate  Flow 


Fig.  3 Semi-log  plots  of  it*  for  adiabatic  wall  boundary  condition 


Backward  Facing  Step 


Fig.  5 Reynolds  stress  profiles  for  the  backward-facing  step  turbulent 
flow  (9:1),  with  data  from  [31] 


Confined  Swirling  Flows 
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Fig.  6 Decay  of  mean  axial  centerline  velocity. 


Confined  Swirling  Flows 


Fig.  8a  Radial  profiles  of  turbulent  intensity(\/i/^) 
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8 - STEP  REACTIONS 


M + O, 

====  o + o 

m + h2 

=====  H + H 

M + H,0 

====  H +0H 

O + H 

====  OH 

H20  + OH 

====  H20  + H 

H2  + OH 

====  H20  + H 

o2  + h 

====  OH  + O 

h2  + o 

====  OH  +H 

k = ATNeXp(-^T) 

A 

N 

E 

0.72000E  + 19 

-1.0000 

117908 

0.55000E+19 

-1.0000 

103298 

0.52000E+22 

-1.5000 

118000 

0.71000E+19 

-1.0000 

0. 

0.58000E+14 

0.0000 

18000 

0.20000E+14 

0.0000 

5166 

0.22000E+15 

0.0000 

16800 

0.75000E+14 

0.0000 

11099 

with  k in  cm3-  mole-1  • s-1  and  E in  cal  • mole-1 


From 


R.C. Rogers  and  Chinitz,  ' Using  a global  hydrogen-air  model 
m turbulent  flow  calcultions  AIAA  J. , vol.  21,  pp.  586-592, 
1983. 
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SSME  Nozzles 


CONTOUR  OF  MACH  NUMBER 


X(U) 


CONTOUR  OF  TEMPERATURE 


Figure  A. 5 Sample  SSME  Nozzle  Flow  Inputs  and  Results  Turbulent, 

8-Step  Kinetics. 

ISP  * 452.78  sec.  Exp.  ISP=453.3sec 
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6'Ai 


distance  from  sub— domain  inlet  aions  the  wall,  inches 


Fig.  17(a),  Pressure  levels  along  the  wall  near  the  nozzle  exit  using  the 
ASM  and  k - e models 


distance  from  sub— domain  inlet  along  the  wall,  inches 


Fig.  17(b),  Effects  of  wall  temperature  on  the  wall  pressure 
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Summaries 
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3-D  Calculations  Are  Desirable 


Summaries 


3-D  Calculations  Are  Desirable 


N95-  23443  ^3/ -jji 


ADVANCED  TURBULENCE  MODELS  FOR  TURBOMACHINERY 


Ali  H.  Hadid,  Michele  E.  DeCroix,  and  Munir  M.  Sindir 
Rocketdyne  Division,  Rockwell  International 


93 
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ABSTRACT 

Development  and  assessment  of  the  single-time-scale  k-e  turbulence  model 
with  different  near-wall  treatments  and  the  multi-scale  k-e  turbulence  model  for 
rotating  flows  are  presented.  These  turbulence  models  are  coded  as  self- 
contained  module  decks  that  can  be  interfaced  with  a number  of  CFD  main  flow 
solvers.  For  each  model,  a stand-alone  module  deck  with  its  own  formulation, 
discretization  scheme,  solver  and  boundary  condition  implementations  is 
presented.  These  satellite  decks  will  take  as  input  (from  a main  flow  solver)  the 
velocity  field,  grid,  boundary  condition  specifications  and  will  deliver  turbulent 
quantities  as  output.  These  modules  were  tested  as  a separate  entities  and 
although  many  logical  and  programming  problems  were  overcome  only  wider 
use  and  further  testing  can  render  the  modules  sufficiently  "fool  proof". 
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DEVELOPMENT  OF  A MODULAR  FORMAT  FOR 
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April  20-22, 1993 

NASA  Marshall  Space  Flight  Center 


TURBULENCE  MODELS  TO  BE  ASSESSED 


s« 

CM* 


CO 

-I 

LU 

o 

CO  g 
Q 2 

7^  CO 
O CO  ^ 
Z HI  5 
DC  (o 
UJhT 
DC  CO  £ 


CO 

_l 

LU 

O 

oo 


<2 

£co 
m co  ^ 

HI  LU  5 

< co  S 


LU  jfJ 

PS 

Q Q- 


< o 


T HI 
Q DC 
O < 


LU 

* Q 

o g 

hi 

Q ^ 


LU  * 

-l  o 


ZD  LU 
G G 


is 


< y 
o * 


_i  < 

o o 


5 W 
co  f= 


— ^ 

P 2 

HI  O 


s f— 

>■  LU 


co  2 

x 2 
< > 


I-  CM 


a 


• NEAR-WALL  TREATMENTS  WILL  INCLUDE  (WHERE  APPROPRIATE)  WALL  FUNCTIONS, 
MULTILAYER  MODELS,  AND  LOW-REYNOLDS  NUMBER  APPROXIMATIONS 


TURBULENCE  MODEL  DECK  STRUCTURE  AND 
INTEGRATION  WITH  NAVIER-STOKES  SOLVER 

MODULES  ARE  BASED  ON  THE  PHENOMENOLOGICAL  SINGLE  POINT  TURBULENCE 
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1 | USER  PROVIDED 

Rockwell  International 

Rocketdyne  Division  KSSSNNJ  ROCKETDYNE  PROVIDED 


SINGLE-SCALE  k-e  MODEL  (CONT'D) 
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Rockwell  International 

Rocketdyne  Division 
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STRONGLY  IMPLICIT  SOLVER 


KEMOD-1  MODULE  DECK  (CONT’D) 
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PARTICULAR  PROBLEM 


KEMOD-1  MODULE  DECK  (CONT'D) 

SUBROUTINE  TWOLAY 

CALLED  IF  THE  TWO-LAYER  OR  THE  LOW-REYNOLDS  NUMBER 
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MOMENTUM  EQUATION  SOLVER  OF  THE  MAIN  ROUTINE.  IT  UPDATES 
THE  FLUX  SOURCE  TERM  OF  THE  DISCRETIZED  MOMENTUM  EQUATION 
DUE  TO  WALL  SHEAR  STRESSES 


KEMOD-1  MODULE  DECK  (CONT’D) 


START 


READ  GEOMETRY  DATA 
GRID  NODES 


CALL  GRIDG  CALCULATE 
Fx,  Fy,  ARE,  VOL 


CALL  CALCE  (TE,  1) 
(TURB  KINETIC  ENERGY  K) 


MODPHI 
FOR  IDIR  - i 

SOLSIP 


IF(LAY2  • OR  »LRE) 


TWO  LAY 


CALL  CALCE  (ED,  2) 
ENERGY  DISSIPATION 
EQUATION  (c) 


MODVIS 


7 

MMDPHI  ] 
FOR  IDIR  = 2| 



SOLSIP 

KEMOD  FLOW  CHART 


Rockwell  International 

Rocketdyne  Division 


CFD  93  0 1 3-01 2/D  t/AHH 
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KEMOD-1  MODULE  DECK  (CONT’D) 


y START 

! 

GEOMETRY  AND  INITIAL 
CONDITIONS 


ITER  = 0 

ITERATICN  COUNTER 


ITER  = ITER  t 1 


U,  V MAIN 
FLOW  SOLVER 


CALL  OTHER  ROUTINES 
FOR  W-EOUATION, 
PRESSURE-CORRECTION 
OR  OTHERS 


AFTER  ASSEMBLING 
COEFFICIENTS  AND 
BEFORE  CALLING 
SOLVER  FOR  U AND  V 

CALL  MCDIF' 


CALL  KEMCD  TO 
CALCULATE 
kANDe 


CHECK  MAXIMUM 
RESIDUAL 

ERROR  OF  EQUATIONS 


NQ  IF  ITER  *GE  • MAXIMUM  ITERATION 
^ OR  RESIDUAL  ERROR  IS 

SMALLER  THAN  PRESET  VALUE 


OUTPUT 


KEMOD-1  INTERFACE  WITH  A MAIN  SOLVER 


Rockwell  International 

Rocketdyne  Division 


CF0  93  013-01M31/AHH 
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MULTI-SCALE  k-e  MODEL  (CONT'D) 
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Rockwell  International 


KEMOD-2  MODULE  DECK 


SUBROUTINE  CALCE  ASSEMBLES  THE  COEFFICIENTS  AND  SOURCE 
TERMS  FOR  THE  DISCRETIZED  Kp,  ep,  Kj,  TRANSPORT  EQUATIONS 


KEMOD-2  MODULE  DECK  (CONT’D) 


r 


KEMOD-2  FLOW  CHART 


Rockwell  International 

Rocketdyne  Division 


CFD  93  0 1 3-01 7/0 1/AHH 
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KEMOD-2  MODULE  DECK  (CONT'D) 


GEOMETRY  AND  INrTIAL 
CONDmONS 


ITER  * 0 

ITERATION  COUNTER 


ITER  a ITER  + 1 


Uf  V MAIN 
FLOW  SOLVER 


CALL  OTHER  ROUTINES 
FOR  W-EGUAT10N, 
PRESSURE-CORRECTION 
OROTHERS 


pAl^ULAlgl 


CHECK  MAXIMUM 
RESIDUAL 

ERROR  OF  EQUATIONS 


NO  IF  ITER  -GE  • MAXIMUM  ITERATION 
^ OR  RESIDUAL  ERROR  IS 

SMALLER  THAN  PRESET  VALUE 


OUTPUT 


KEMOD-2  INTERFACE  WITH  MAIN  SOLVER 


Rockwell  International 

Rocketdyne  Division 


CF0  93  013-018/D1/AHH 
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ROBACK  AND  JOHNSON  RESULTS 
TANGENTIAL  VELOCITY  PREDICTIONS  AT  X = 0.025  M 
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SUMMARY 
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WORK  ON  3-D  MODULES  TO  BEGIN  AS  SCHEDULED  (FY  '94) 


N95- 23444 


/99a-  1 1 ?OcD  V 

NUMERICAL  COMPUTATION  OF  AERODYNAMICS  AND  HEAT  TRANSFER  IN  A 
TURBINE  CASCADE  AND  A TURN-AROUND  DUCT  USING  ADVANCED 

TURBULENCE  MODELS 

B.  LAKSHMINARAYANA  AND  J.  LUO 

The  Pennsylvania  State  University 
Department  of  Aerospace  Engineering 
University  Park,  PA  16802 

The  objective  of  this  research  is  to  develop  turbulence  models  to  predict  the  flow  and  heat 
transfer  fields  dominated  by  the  curvature  effect  such  as  those  encountered  in  turbine  cascades  and 
turn-around  ducts. 

A Navier-Stokes  code  has  been  developed  using  an  explicit  Runge-Kutta  method  with  a two 
layer  k-e/ARSM  (Algebraic  Reynolds  Stress  Model),  Chien’s  Low  Reynolds  Number  (LRN)  k-e 
model  and  Coakley’s  LRN  q-  model.  The  near  wall  pressure  strain  correlation  term  was  included 
in  the  ARSM.  The  formulation  is  applied  to  Favre-averaged  N-S  equations  and  no  thin-layer 
approximations  are  made  in  either  the  mean  flow  or  turbulence  transport  equations.  Anisotropic 
scaling  of  artificial  dissipation  terms  was  used.  Locally  variable  timestep  was  also  used  to  improve 
convergence.  Detailed  comparisons  were  made  between  computations  and  data  measured  in  a turbine 
cascade  by  Arts  et  al.  at  Von  Karman  Institute.  The  surface  pressure  distributions  and  wake  profiles 
were  predicted  well  by  all  the  models.  The  blade  heat  transfer  is  predicted  well  by  k-e/ARSM 
model,  as  well  as-  the  k-  e model.  It’s  found  that  the  onset  of  boundary  layer  transition  on  both 
surfaces  is  highly  dependent  upon  the  level  of  local  freestream  turbulence  intensity,  which  is  strongly 
influenced  by  the  streamline  curvature. 

Detailed  computation  of  the  flow  in  the  turn  around  duct  has  been  carried  out  and  validated 
against  the  data  by  Monson  as  well  as  Sandborn.  The  computed  results  at  various  streamwise  locations 
both  on  the  concave  and  convex  sides  are  compared  with  flow  and  turbulence  data  including  the 
separation  zone  on  the  inner  well.  The  k-e/ARSM  model  yielded  relatively  better  results  that  the 
two-equation  turbulence  models.  A detailed  assessment  of  the  turbulence  models  has  been  made  with 
regard  to  their  applicability  to  curved  flows. 


?' 
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NUMERICAL  COMPUTATION  OF  AERODYNAMICS  AND 
HEAT  TRANSFER  IN  A TURBINE  CASCADE  AND  A 
TURN-AROUND  DUCT  USING  ADVANCED 
TURBULENCE  MODELS* 

B.  Lakshminarayana  and  J.  Luo 

The  Pennsylvania  State  University 
Department  of  Aerospace  Engineering 
University  Park,  Pennsylvania 

OBJECTIVE: 


TO  DEVELOP  TURBULENCE  MODELS  TO  PREDICT 
FLOW  AND  HEAT  TRANSFER  FIELDS  IN 
TURBOMACHINERY  INCLUDING  CURVATURE, 
ROTATION  AND  HIGH  TEMPERATURE  EFFECTS 


OUTLINE: 

• INTRODUCTION 

• NUMERICAL  TECHNIQUE 

• TURBULENCE  MODELS 

• FLOW  AND  HEAT  TRANSFER  FIELD  IN  A HIGH 
MACH  NUMBER  TRANSONIC  TURBINE  CASCADE 

• FLOW  FIELD  IN  A TURN-AROUND  DUCT 

• CONCLUSIONS 


•SPONSORED  BY  NASA  HUNTSVILLE  WITH  USA  GRIFFIN  AS  THE  TECHNICAL 
MONITOR 
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Relative  velocities,  constant  rotation  rate  about  X axis,  co.  Averaged  quantities. 

ai  o)V 

Energy,  ®°’£  2 ' 2 .rothalpy  constant  along  streamlines  for  inviscid  steady  state. 


TECHNIQUES 


1.  RK2D  code  : 

* 2-D  Navier-Stokes  code, 

Conservative,  compressible  formulation 

* Favre- Averaged  Mean  and  Turbulence 
equations 

* 4-stage  explicit  Runge-Kutta  scheme 

* 2nd  and  4th  order  artificial  dissipation 
(with  eigenvalue  and  local  velocity  scaling) 

* Coupled  with  compressible  Low- 

Reynolds  number  K-e  model , q-co  model , 
ARSM  , NLSM  ( Nonlinear  -stress  model) , 
AHFM(Algebraic  Heat  Flux  model) 

* Characteristic  boundary  conditions,  H 
grids  (generated  by  a combined  algebraic 
and  elliptic  method  to  keep  smoothness  and 
orthogonality  near  the  wall) 
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2.  TEXSTAN  code 

* 2-D  boundary  layer  code  developed  by 
Crawford 

* Extension  of  STAN  5,  Patankar- 
Spalding  numerical  scheme 

* Include  7 differential  two-equation 
turbulence  models  (Jones-Launder, 
Chien,Lam-  Bremhorst  , etc.)  and  mixing 
length  model 
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RSM  ( Reynolds  Stress  Model 


(Gibson  & Launder  1978) 
Reynolds  stresss  transport  equation  : 


5U;U; 


ux  sk1= -wj>k  - Wi,k +>i,j + v 


5 Pui?  piij _ du.u]  5u.  3uj 

^ — uiuiu],+ — ^5..  +i— i5.,  -y— . 1 1 _2v l I 

^xk  J k P k P Jk  dxk  dxh  d\. 


i.e.,  Cy-Dy-Pj. 

where  Ey=|e  (Dissipation) 

‘t’ij  =<fiji+*t*ij2+V,w+<t'ij2,w  (Pressure-strain 
correlation) 

‘t’iji =-cif  (5i“j_fkV  (Retum-to-isotropy 
part) 

V =-c2(pij  -|pk8ij)  (Rapid  part) 
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<,ijw4=cl|(ukuninknm5ij-|uiuknknj-|uiuknkni)fn 

(Near-wall  term) 

V,2  = C2((l)km,2nknm8ij  “ 2<t>jMnkni)f» 

(Near-wall  term) 

rn=k3/2/(2.5SxnE)  (x„  is  the  distance  normal 
to  the  wall) 

Constants  : cx=1.8,  c2=0.6,  c;  =0.5,  cj 
=0.3 
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ARSM  (Algebraic  Reynolds  Stress  Model) 
ARSM  assumption: 


U-U  ■ U;U; 

cij _ Du = -ir(Ck  - Dk} = -ir(pk  -£) 


u.u 


U-U. 


=>  Cij-Dij  = -irl(Ck-Dk)  = ^i(Pk-e) 


U:U: 


U-U- 


=>  cij-Dij=-Tr(Cx-Dk)=-ir(Pk-£) 


where 


P.  =— u u U.  . 

k ij  i,j 


pij=-uiukuj>k-*vI^ui,k 


<t>ijl;W  and  <t>ii2.w  as  in  RSM 


ij2,w 
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NLSM  (Nonlinear  Stress  Model 
(Shih  , Zhu  & Lumley  1992) 


Reynolds  stress  : 


uiuj  = 3k5ij"Vt(UiJ  + Uj.i) 

+ A^VS(U^UM  + UMUK,i-§7t5ij) 


c„  k> 
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where 


5t  = Ui,jUj,i 

^ = uijui,j 

fc2 

v.  = chV 


^ Aj+ri  + a^ 


$ = *£} 
^ e 


Q = (2Q*ja*j)2 

Q*  = (U.  .-U-.)/2 
u v i»j  j>i' 

r|  = — S 

1 e 


S = (2S..S..)2 

v ijij; 


S..=(U.  . + U..)/2 
ij  v i,j  j,i' 

Constants: 


CT1 

Ct2 

C,2 

A! 

a 

-4 

13 

-2 

1.25 

0.9 

1000 
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turbulence  Intensity 


Fig.  Validation  of  the  ARSM  model: 
Turbulence  intensity  profile  in  the  flat-plate 
turbulent  boundary  layer  : experiment  by 

Klebanoff;  computation  by  ARSM 
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1 80-riegree  TURN  AROUND  DUCT (TAD 
Geometry  & Grid  : 


I 


EARLIER  RESEARCH  ON  TAD  FT  ,OW 


• Measurements: 

*Sandbom  (1988),  Sandbom  and  Shin 
(1989)  (Water  flow,  Re=7xio4-  s*io5  (Re 
based  on  duct  height  and  bulk  velocity) 

*Monson,  Seegmiller  , McConnaughey 
& Chen  (1989,1990)  (Air  flow,  Re=io5,io6) 

*Sharma  etal  (1987)  (Axisymmetric 
TAD  air  flow,  Re=io5) 
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Earlier  Computations: 


*Chen  and  Sandbom  (1986)  (K-e  and 
curvature-corrected  K-e  ) 

*Monson,  Seegmiller  & McConnaughey 

(1989,1990)  (mixing-length  & K-e  with 
curvature-correction) 

*Avva  etal  (1990)  (High  Re  and  Low- 
Re  K-e ) 

*Gallardo  & Lakshminarayana  (1993) 
(curvature -modified  K-e ) 

• Agreement  in  above  computations  are 
not  satisfactory. 
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A data,  inner  wall  (IW) 
O data  .outer  wall  (OW) 

K-E  (IW) 

- - • K-E  (OW) 

■ NLSM  (IW) 

NLSM  (OW) 

ARSM  (IW) 

ARSM  (OW) 


Fig.l  Static  pressure  coefficient  on  turnaround 


duct  inner  and  outer  walls 
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Fig.  2 Skin  friction  coefficient  on  turnaround 
duct  inner  and  outer  walls 
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Fig.3  (a)  Longitudinal  velocity  in  turnaround 
duct,  x/H=-4 
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Fig.3(b)  Turbulent  kinetic  energy  profile 
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Fig.5(b)  Longitudinal  velocity  in 
turnaround  duct,  theta=90  deg. 
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Fig.5(c)  Turbulent  kinetic  energy 
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Fig.5(d)  Turbulent  shear  stress 
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Fig.5  (e)  Comparison  of  K and 
0.75*(uu+vv) 
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Fig. 7(a)  Longitudinal  velocity  in  turnaround 
duct,  x/H  = 2 (downstream  of  turn) 
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Fig.  7(c)  Turbulent  shear  stress 
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WT  Turbine  Nozzle  Guide  Vane  Cascade 


* Measurement  by  Arts  etal  (1990)  at 
Von  Karman  Institute 

* M(inlet)=0.15,  M(outlet)=0.7  to 
1.11,  Re=0.5  - 2 x io«  , To=415k,  T(wall)  = 

300k 


* Geometry  and  grid 


3 . Smoothing:  Typical  values  for 

2nd  and  4th  order  dissipation 
taken  as  2% -3%  and  3% -4%  for  the 
turbine  cascade  computations. 

4 . Table:  Computed  cases  : 


cases 

Mur228 

Mur224 

Mur239 

Pol(bar) 

0.915 

0.909 

3.387 

T01  (K) 

403 

403 

412 

Re, 2 

1 0.6E  + 6 

0.6E+6 

[2.1  E+6 

Tujn(%) 

1 

6 

6 

T(wall)  = 300k  for  all  these  cases. 
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Po1-Po2  (mmHg) 


0 20  40  60  80 

S (mm) 


Fig.10  Computed  and  measured  wakes  for 
mur043 
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(w/m2k) 


(W/m2k) 
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CONCLUSIONS 


*A  two-dimensional  Navier-Stokes 
code  has  been  developed  using  an  explicit 
Runge-Kutta  method  incorporating  the 
ARSM  model,  NLSM  model , Chien’s  LRN 

k-e  model  and  Coakley's  LRN  q-co  model. 

*The  surface  pressure  distributions  and 
wake  profiles  of  a transonic  turbine  cascade 

were  predicted  well  by  the  k-e  , ARSM 

and  q-G)  models.  The  heat  transfer  on 
suction  surfaces  were  predicted  well  by  the 

k-e  and  ARSM  models. 

*The  heat  transfer  on  pressure  side  for 
one  case  (MTJR239)  was  underpredicted. 
This  was  caused  by  the  underprediction  of 
mainstream  turbulence  level  , which 
strongly  influences  the  transition  location. 

*The  boundary  layer  code  predicts  the 
heat  transfer  on  pressure  surface  well,  but  it 
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does  not  capture  the  transition  on  suction 
surfaces  for  all  the  3 cases. 

*The  wall  damping  function  (r^)  in 
Chien's  model  was  modified  to  yield 
improved  prediction  for  flow  under  adverse 
pressure  gradient. 

*For  TAD  flow,  good  predictions  have 
been  obtained  for  the  surface  pressure 
distribution  and  skin  friction  coefficients. 
The  ARSM  model  yields  better  prediction 

than  NLSM  and  k-8  models,  for  both  the 
mean  and  the  turbulence  quantities. 

*The  near  wall  "echo"  term  <t>ijw  is  not 

correctly  modelled  for  strongly  curved 
flows,  especially  near  the  concave  surface. 

*For  more  accurate  prediction  of 
strongly  curved  TAD  flows,  it  may  be 
necessary  to  use  existing  RSM,  including 
the  modeling  for  the  wall  region. 
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Liquid  Propellant  Rocket  Engine  Combustion  Simulation 
with  a Time-Accurate  CFD  Method 


Y.S.  Chen,  H.M.  Shang  and  Paul  Liaw 
Engineering  Sciences,  Inc. 

Huntsville,  Alabama 

and 
J.  Hutt 

NASA  George  C.  Marshall  Space  Flight  Center 
NASA/MSFC,  Alabama 

APSTRAC1 

Time-accurate  computational  fluid  dynamics  (CFD)  algorithms  are 
among  the  basic  requirements  as  an  engineering  or  research  tool  for  realistic 
simulations  of  transient  combustion  phenomena,  such  as  combustion 
instability,  transient  start-up,  etc.,  inside  the  rocket  engine  combustion 
chamber.  A time-accurate  pressure  based  method  is  employed  in  the  FDNS 
code  for  combustion  model  development.  This  is  in  connection  with  other 
program  development  activities  such  as  spray  combustion  model  development 
and  efficient  finite-rate  chemistry  solution  method  implementation.  In  the 
present  study,  a second-order  time-accurate  time-marching  scheme  is 
employed.  For  better  spatial  resolutions  near  discontinuities  (e.g.  shocks, 
contact  discontinuities),  a 3rd-order  accurate  TVD  scheme  for  modeling  the 
convection  terms  is  implemented  in  the  FDNS  code.  Necessary  modification 
to  the  predictor/multi-corrector  solution  algorithm  in  order  to  maintain  time- 
accurate  wave  propagation  is  also  investigated.  Benchmark  1-D  and  multi- 
dimensional test  cases,  which  include  the  classical  shock  tube  wave 
propagation  problems,  resonant  pipe  test  case,  unsteady  flow  development  of 
a blast  tube  test  case,  and  H2/02  rocket  engine  chamber  combustion  start-up 
transient  simulation,  etc.,  are  investigated  to  validate  and  demonstrate  the 
accuracy  and  robustness  of  the  present  numerical  scheme  and  solution 
algorithm. 
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MODELS  FOR  STEADY-STATE  AND  TRANSIENT 
FLOW  APPLICATIONS  IS  THE  FINAL  GOAL 


BASIC  BUILDING  BLOCKS 
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GAS  PHASE  GOVERNING  EQUATIONS 
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VALIDATION  CASES 
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Convergence  Acceleration  of  Implicit  Schemes  ^ ^ ' / / 9 006, 
in  the  Presence  of  High  Aspect  Ratio  Grid  Cells 


P.  E.  0.  Buelow,  S.  Venkateswaran  and  C.  L.  Merkle 
Propulsion  Engineering  Research  Center 
The  Department  of  Mechanical  Engineering 
The  Pennsylvania  State  University 
University  Park,  PA  16802. 


The  performance  of  Navier-Stokes  codes  are  influenced  by  several  phenomena.  For 
example,  the  robustness  of  the  code  may  be  compromised  by  the  lack  of  grid  resolution, 
by  a need  for  more  precise  initial  conditions  or  because  all  or  part  of  the  flowfield 
lies  outside  the  flow  regime  in  which  the  algorithm  converges  efficiently.  A primary 
example  of  the  latter  effect  is  the  presence  of  extended  low  Mach  number  and/or  low 
Reynolds  number  regions  which  cause  convergence  deterioration  of  time  marching 
algorithms.  Recent  research  into  this  problem  by  several  workers  including  the  present 
authors  has  largely  negated  this  difficulty  through  the  introduction  of  time-derivative 
preconditioning.  In  the  present  paper,  we  employ  the  preconditioned  algorithm  to 
address  convergence  difficulties  arising  from  sensitivity  to  grid  stretching  and  high  aspect 
ratio  grid  cells. 

Strong  grid  stretching  is  particularly  characteristic  of  turbulent  flow  calculations 
where  the  grid  must  be  refined  very  tightly  in  the  dimension  normal  to  the  wall,  without 
a similar  refinement  in  the  tangential  direction.  High  aspect  ratio  grid  cells  also  arise  in 
problems  that  involve  high  aspect  ratio  domains  such  as  combustor  coolant  channels. 

In  both  situations,  the  high  aspect  ratio  cells  can  lead  to  extreme  deterioration  in 
convergence.  It  is  the  purpose  of  the  present  paper  to  address  the  reasons  for  this 
adverse  response  to  grid  stretching  and  to  suggest  methods  for  enhancing  convergence 
under  such  circumstances. 

Numerical  algorithms  typically  possess  a maximum  allowable  or  optimum  value 
for  the  time  step  size,  expressed  in  non-dimensional  terms  as  a CFL  number  or  von- 
Neumann  number  (VNN).  In  the  presence  of  high  aspect  ratio  cells,  the  smallest 
dimension  of  the  grid  cell  controls  the  time  step  size  causing  it  to  be  extremely  small, 
which  in  turn  results  in  the  deterioration  of  convergence  behaviour.  For  explicit  schemes, 
this  time  step  limitation  cannot  be  exceeded  without  violating  stability  restrictions  of 
the  scheme.  Oh  the  other  hand,  for  implicit  schemes,  which  are  typically  unconditionally 
stable,  there  appears  to  be  room  for  improvement  through  careful  tailoring  of  the  time- 
step  definition  based  on  results  of  linear  stability  analyses.  In  the  present  paper,  we 
focus  on  the  central-differenced  alternating  direction  implicit  (ADI)  scheme.  The 
understanding  garnered  from  this  analyses  can  then  be  applied  to  other  implicit 
schemes. 

In  order  to  systematically  study  the  effects  of  aspect  ratio  and  the  methods  of 
mitigating  the  associated  problems,  we  use  a two  pronged  approach.  We  use  stability 
analyses  as  a tool  for  predicting  numerical  convergence  behavior  and  numerical 
experiments  on  simple  model  problems  to  verify  predicted  trends.  Based  on  these 
analyses,  we  determine  that  efficient  convergence  may  be  obtained  at  all  aspect  ratios 
by  getting  a combination  of  things  right.  Primary  among  these  axe  the  proper  definition 
of  the  time  step  size,  proper  selection  of  viscous  preconditioner  and  the  precise  treatment 
of  boundary  conditions.  These  algorithmic  improvements  are  then  applied  to  a variety  of 
test  cases  to  demonstrate  uniform  convergence  at  all  aspect  ratios. 
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Philosophy  of  Grid  Aspect  Ratio  Study 


• Assessment  of  High  Aspect  Ratio  Problem 

— Disparate  propagation  speeds  in  X and  Y 

• Stability  Theory 

— Scalar  Convection-Diffusion  Equation 

— Euler  Equations 

— Navier-Stokes  Equations 

• Numerical  Convergence  Studies 

— Simple  Model  Problems 

— Realistic  Flow  Problems 

• Improved  Algorithm  to  Provide  Aspect  Ratio  Control 

— Precise  Time-Step  Definition 

— Viscous  Preconditioning 

— Boundary  Condition  Implementation 
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The  Navier-Stokes  Equations 


r$Q  v 
dt 


dE  dF_ 
+ dx  dy 


= H + L(QV) 


• Solution  Vector 


Qv  = ( P,  u,  v , T )T 
• Preconditioning  Matrix 


1 fee2 

0 

0 

0 \ 

u/ec2 

p 

0 

0 

v fee1 

0 

p 

0 

\W-W) 

ec2  ~~  1 

pu 

pv 

pej 

• Parameter  e 

— Activates  Inviscid  and  Viscous  Preconditioning 

— Value  Depends  on  Local  Mach  Number  and  Cell 
Reynolds  Number 
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Numerical  Solution  Procedure 


Central-Differenced  ADI  Algorithm 


_ dA  d n 
S + -r —Rxx 

OX  ox 


r-1 


' dB  d „ ■ 

+ dy  dy^. 


A<?„  = -TC" 


• Approximate  factorization  errors  control  convergence 
behavior. 

• Optimum  CFLu+c  is  typically  between  1 and  10. 

— Other  inviscid  and  viscous  time  scales  are  opti- 
mized by  the  preconditioning  matrix. 
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Time-Step  Definition 


• Local  Time-Stepping  or  Constant  CFL  Condition 

Max  ( CFLX , CFLy  ) = CFL 

• For  high  aspect  ratios,  CFLX  and  CFLy  become  dis- 
parate 

CFLy  = CFL , CFLX  = CFL  JAR 
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Euler  Stability  Analysis 

• Aspect  Ratio  (AR)  of  100 

CFLy  = 1 
CFLX  = 0.01 

l/)y 

0 

CFLy  = 10 
CFLX  =0.1 

U) 

o fy 
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New  Time-Step  Definition 


Conclusions  from  Stability  Analysis: 

— Min-CFL  Preferable  to  Ma x-CFL 

— Efficient  Convergence  at  all  AR 

Minimum-CFL  Definition 


Min  ( CFLX , CFLy  ) = CFL 


For  high  aspect  ratios, 


CFLy  = CFL  * AR, 


CFLX  = CFL 


Implementation  of  Boundary  Conditions 


• Extrapolation  vs  Characteristic 

— Both  work  well  for  small  CFV s 

— Characteristics  usually  superior  at  high  CFV s 

• Proper  MOC  Implementation: 

— Implicit  procedures 

— Boundary  conditions  applied  before  approximate  fac- 
torization 

— Consistent  order  of  accuracy:  LHS  / RHS 
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No.  of  Iterations  to  Machine  Accuracy 


High  Aspect  Ratio  Convergence — Inviscid  Duct 


le+05  f 


le+04  \r 


1000  t- 


100  L 
0.1 


• • | • i * I "■  i 1 r ‘t  | i i i 


1 10  100  1000 
Grid  Aspect  Ratio 
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Navier-Stokes  Analysis 


• Parameter  e controls  low  Re  number  convergence 

— Viscous  terms  limit  time  step  at  high  AR 

— 6 chosen  to  optimize  inviscid  and  viscous  modes  simul- 
taneously 

• Obvious  choice:  e = / (Ma x-CFL,  Ma x-VNN) 

• Scalar  Stability  Results:  e = / (Min-CFL,  Min -VNN) 
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Navier-Stokes  Analysis 


• Min-CFL,  Min-F NN  Stability  Result 
AR  = 1000 


VNNX  = 1,  VNNy  = 1 x 106 
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Navier-Stokes  Analysis 


• Conclusions  from  Stability  Results: 

— Vector  system  different  from  scalar  equation 

— Approximate  factorization  error  CFLX  *VNNy  limits 
convergence 

• Viscous  preconditioner,  e = / (Min-CFL,Max-ViViV) 

— Maintains  Min-CFL  for  ‘inviscid’  modes 

— Uses  traditional  Max-ViVJV  definition  for  ‘viscous’  modes 
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Navier-Stokes  Analysis 

• Min-CFL,  Ma x-VNN  Stability  Result 
AR  = 1000 


CFLX  = 1,  CFLy  = 1000 
VNNX  = 1 x 10"6,  VNNy  = 1 
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No.  of  Iterations  to  Machine  Accuracy 


High  Aspect  Ratio  Convergence — Viscous  Duct 


Grid  Aspect  Ratio 
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RADIAL  COORDINATE 


Hydrogen/Oxygen  Shear-Layer 
Stretched  Grid  and  Solution 
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Hydrogen/ Oxygen  Shear-Layer 
Convergence 


No.  of  Iterations 
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Norm  of  the  Residual 


High  Reynolds  Number  Boundary  Layer 
Convergence — Stretched  Grid 
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High  Reynolds  Number  Boundary  Layer 

Blasius  Solution 


u/U 
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High  Reynolds  Number  Boundary  Layer 
Turbulent  Velocity  Profile 
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Turbulent  Nozzle  Computation 
Stretched  Grid  and  Solution 
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L2  Norm  of  the  Residual 


Turbulent  Nozzle  Computation 
Convergence 


No.  of  Iterations 
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Heat  Flux,  W/ 


Turbulent  Nozzle  Computation 
Wall  Heat  Flux 
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Conclusions — High  Aspect  Ratio  Study 


• High  Aspect  Ratio  Analysis 

— Stability  Theory 

— Numerical  Convergence  Studies 

• Convergence  Control: 

— Min -CFL  Time  Step 

— Ma x-VNN  Viscous  Preconditioner 

— Correct  implementation  of  boundary  conditions 

• Uniform  Convergence  Demonstrated  for  All  AR’s 

— Above  issues  addressed  in  combination 

— Efficient  convergence  for  variety  of  test  cases 
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Conclusions  (ContcL) 


• Present  results  are  for  two-dimensional  central-differenced 
ADI  scheme 

• Explicit  Schemes: 

— Optimum  time  step  causes  poor  convergence  at  high 
Ai R’s 

• Upwind  Schemes  Also  Suffer  at  High  AR’s 

— Present  improvements  may  be  incorporated 

— Rich  variety  of  approximate  factorization  methods 

• Three-dimensional  computations: 

— ADI  scheme  is  conditionally  stable 

— Two  kinds  of  high  aspect  ratio  grids 

— Algorithmic  improvements  appear  promising 
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A TIME- ACCURATE  FINITE  VOLUME  METHOD 
ALL  FLOW  VELOCITIES 

S.-W.  Kim 

Resident  Research  Associate 
CFD  for  Space  Branch, 

NASA  Lewis  Research  Center,  MS  5-11, 

Cleveland,  Ohio  44135 

A finite  volume  method  to  solve  the  Navier-Stokes  equations  at  all  flow  velocities  (e.g., 
incompressible,  subsonic,  transonic,  supersonic  and  hypersonic  flows)  is  presented.  The 
numerical  method  is  based  on  a finite  volume  method  that  incorporates  a pressure-staggered 
mesh  and  an  incremental  pressure  equation  for  the  conservation  of  mass.  Comparison  of  three 
generally  accepted  time-advancing  schemes,  i.e..  Simplified  Marker-and-Cell  (SMAC), 
Pressure-Implicit- Splitting  of  Operators  (PISO),  and  Iterative-Time- Advancing  (TTA)  scheme, 
are  made  by  solving  a lid-driven  polar  cavity  flow  and  self-sustained  oscillatory  flows  over 
circular  and  square  cylinders.  Calculated  results  show  that  the  ITA  is  the  most  stable 
numerically  and  yields  the  most  accurate  results.  The  SMAC  is  the  most  efficient 
computationally  and  is  as  stable  as  the  ITA.  It  is  shown  that  the  PISO  is  the  most  weakly 
convergent  and  it  exhibits  an  undesirable  strong  dependence  on  the  time-step  size.  The 
degenerated  numerical  results  obtained  using  the  PISO  is  attributed  to  its  second  corrector  step 
that  cause  the  numerical  results  to  deviate  further  from  a divergence  free  velocity  field.  The 
accurate  numerical  results  obtained  using  the  ITA  is  attributed  to  its  capability  to  resolve  the 
nonlinearity  of  the  Navier-Stokes  equations. 

The  present  numerical  method  that  incorporates  the  ITA  is  used  to  solve  an  unsteady 
transitional  flow  over  an  oscillating  airfoil  and  a chemically  reacting  flow  of  hydrogen  in  a 
vitiated  supersonic  airstream.  The  turbulence  fields  in  these  flow  cases  are  described  using 
multiple-time-scale  turbulence  equations. 

For  the  unsteady  transitional  over  an  oscillating  airfoil,  the  fluid  flow  is  described  using 
ensemble-averaged  Navier-Stokes  equations  defined  on  the  Lagrangian-Eulerian  coordinates. 
It  is  shown  that  the  numerical  method  successfully  predicts  the  large  dynamic  stall  vortex 
(DSV)  and  the  trailing  edge  vortex  (TEV)  that  are  periodically  generated  by  the  oscillating 
airfoil.  The  calculated  streaklines  are  in  very  good  comparison  with  the  experimentally 
obtained  smoke  picture.  The  calculated  turbulent  viscosity  contours  show  that  the  transition 
from  laminar  to  turbulent  state  and  the  relam inarization  occur  widely  in  space  as  well  as  in 
time.  The  ensemble-averaged  velocity  profiles  are  also  in  good  agreement  with  the  measured 
data  and  the  good  comparison  indicates  that  the  numerical  method  as  well  as  the  multiple- 
time-scale  turbulence  equations  successfully  predict  the  unsteady  transitional  turbulence  field. 

The  chemical  reactions  for  the  hydrogen  in  the  vitiated  supersonic  airstream  are  described 
using  9 chemical  species  and  48  reaction-steps.  Consider  that  a fast  chemistry  can  not  be  used 
to  describe  the  fine  details  (such  as  the  instability)  of  chemically  reacting  flows  while  a 
reduced  chemical  kinetics  can  not  be  used  confidently  due  to  the  uncertainty  contained  in  the 
reaction  mechanisms.  However,  the  use  of  a detailed  finite  rate  chemistry  may  make  it  difficult 
to  obtain  a fully  converged  solution  due  to  the  coupling  between  the  large  number  of  flow, 
turbulence,  and  chemical  equations.  The  numerical  results  obtained  in  the  present  study  are  in 
good  agreement  with  the  measured  data.  The  good  comparison  is  attributed  to  the  numerical 
method  that  can  yield  strongly  converged  results  for  the  reacting  flow  and  to  the  use  of  the 
multiple-time-scale  turbulence  equations  that  can  accurately  describe  the  mixing  of  the  fuel 
and  the  oxidant. 
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Multiple-time-scale  turbulence  equations  are  used  in  the  present  unsteady 
transitional  flow  and  the  chemically  reacting  flow  calculations. 
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1.  Iteartive  Time- Advancing  Scheme  (ITA) 
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2.  Simplified  Marker  and  Cell  (SMAC) 
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Solve  eqs.  (1-4)  for  each  time-level.The  velocity  field  strongly  satisfy  the 
conservation  of  mass. 
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3.  Pressure-Implicit  Splitting  of  Operators 
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(pc.+Ar1)  dx, 
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Polar  cavity  flow  starting  from  rest 
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Polar  cavity  flow  starting  from  rest 
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cylinder. 


Strouhal  Numbers  (St  = 2fr/Uo)  for  Flow  over  a Circular  Cylind 
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Ref.  [2] : Braza,  Chassing,  Ha  Minh,  JFM  1986, 
Ref.  [14]:  Eaton,  JFM  1987. 


Calculation  of  Unsteady  Transitional  Flow 
over  Oscillating  Airfoil* 
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* Kim,  Zaman,  and  Panda,  To  be  presented  in  the  ASME  Fluid 
Engineering  Conference,  Washington  D.C.,  June  20-24,  1993. 
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Ensemble  average 
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t/T=0.0  t/T=0.2  t/T=0.4  t/T=0.5  t/T=0.7  t/T=0.9 

(5°u)  (H°u)  (22.5°u)  (25°u)  (19°d)  (7.4°d) 
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Comparison  of  calculated  streaklines  with  smoke  picture. 
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0,  - 
(g)  7.4° 

Turbulent  viscosity  contour  ( A|i/|it  = 10) 
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Chemical  species  concentration  equation 
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Chemical  reactions  for  the  combustion  of  H2  in  a vitiated  supersonic  airstream  are 
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CONCLUSIONS  AND  DISCUSSION 
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The  ITA  can  best  resolve  the  nonlinearity  of  the  Navier-Stokes  Equations. 
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conservation  of  mass,  and 

(iii)  the  Multiple-time-scale  turbulence  equations  that  can  resolve  the 
inequilibrium  turbulence  field. 
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A CONTROLLED  VARIATION  SCHEME  FOR  CONVECTION 
TREATMENT  IN  PRESSURE-BASED  ALGORITHM 


P.  $° 


Wei  Shyy,  Siddharth  Thakur  and  Kevin  I\icker 

Department  of  Aerospace  Engineering,  Mechanics  & Engineering 

Science, 

University  of  Florida 
and 

NASA  Marshall  Space  Flight  Center 


Convection  effect  and  source  terms  are  two  primary  sources  of  difficulties 
in  computing  turbulent  reacting  flows  typically  encountered  in  propulsion 
devices.  The  present  work  intends  to  elucidate  the  individual  as  well  as  the 
collective  roles  of  convection  and  source  terms  in  the  fluid  flow  equations, 
and  to  devise  appropriate  treatments  and  implementations  to  improve  our 
current  capability  of  predicting  such  flows.  A controlled  variation  scheme 
(CVS)  has  been  under  development  in  the  context  of  a pressure-based 
algorithm,  which  has  the  characteristics  of  adaptively  regulating  the 
amount  of  numerical  diffusivity,  relative  to  central  difference  scheme, 
according  to  the  variation  in  local  flow  field.  Both  the  basic  concepts  and 
a pragmatic  assessment  will  be  presented  to  highlight  the  status  of  this 
work. 
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Source  Terms:  Cause  Numerical  Difficulties  Due  to  Different  Length  and 

Time  Scales 


Present  Approach 
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Linear  Steady  Burgers’  Equation 
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Critical  Value  : I P I > 2 
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Model  Problem  II:  Different  Schemes 


1898 


21 


'Pcvs  Proflle 


Solution  profile 
""(not to’scale) 


Model  Problem  II:  Dirichlet  B.C.,  P=100 
Effective  Cell  Peclet  Number 
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Treatment  of  Source  Terms 
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TREATMENT  OF  SOURCE  TERMS:  A LONGITUDINAL 
COMBUSTION  INSTABILITY  PROBLEM 


Besides  convection  terms,  source  terms  (if  present)  in  the  Navier-Stokes 
equations  can  be  strong  enough  to  cause  numerical  difficulties  such  as  a loss 
of  accuracy  in  the  form  of  spurious  oscillations  in  the  solution  profiles  or 
numerical  instability.  This  is  so  because  strong  source  terms  can  be 
sufficiently  stiff  and  the  time  and  length  scales  imposed  by  them  may  not 
be  commensurate  with  those  imposed  by  convection,  for  example.  Thus, 
due  attention  has  to  be  paid  to  the  source  terms  and  not  just  to  the  convection 
terms. 

A one-dimensional  longitudinal  combustion  instability  problem  is  chosen 
which  has  a strong  heat  release  source  term.  The  high  accuracy  TVD  type 
of  convection  treatment  in  a sequential  solver  (second  fig.  clockwise:  top 
right)  is  seen  to  provide  higher  accuracy  than  the  first-order  upwind 
scheme  (first  figure:  top  left),  as  evident  from  the  amplitudes  of  the  ten 
pressure  modeshapes  shown  in  the  viewgraph.  However,  the  TVD  type  of 
convection  treatment  without  any  special  source  term  treatment  yields 
spurious  oscillations  in  modehapes  numbered  5,  6 and  7 (second  figure 
clockwise).  From  the  corresponding  heat  release  modeshapes  (third  fig. 
clockwise),  it  is  clear  that  modes  5, 6 and  7 are  the  modes  of  maximum  heat 
release,  thus  demonstrating  that  when  source  terms  become  stiff  enough 
they  may  lead  to  spurious  oscillations.  This  can  be  resolved  by  increasing 
the  amount  of  numerical  dissipation  in  the  scheme  (by  varying  6)  but  this 
is  accompanied  by  an  overall  smearing  of  solution  profiles.  However, 
special  source  term  treatment  such  as  MacCormack’s  predictor-corrector 
method  or  Strang’s  time-splitting  method  (here,  the  latter)  can  resolve  the 
problem  by  suppressing  any  spurious  oscillations  without  the  need  of  any 
extra  numerical  damping.  This  is  clearly  evident  from  the  bottom  left  plot 
(fourth  fig.  clockwise). 
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Special  Source  Term  Treatment 

Conservation  law  with  a source  term 

+ f(W)x  = 


Treatment 

a)  MacCormack’s  Predictor-Corrector  Method 

b)  Operator  Splitting  (Strang’s  Time-Splitting) 

wn+1  = 5f(A«/2)  Sf(At)  Sq(At/2)  Wn 

where  Sf  represents  the  numerical  solution  operator  for 

wt  + f(w)x  = 0 

and  is  the  numerical  solution  operator  for  the  ODE 


1-D  Combusting  Flow  in  a Duct 

# To  illustrate  the  effect  of  a strong  source  term  (heat 
release)  on  numerical  accurcay 

# To  demonstrate  the  efficacy  of  special  source  term 
treatment  for  a strong  source  term 


3p  + d(?u)  = 0 

dt  dx 


dm  d(mu ) _ dp  1 ^Tll 

dt  dx  dx  Re^  dx 


dE  + d(Eu) 
dt  dx 


d(pu) 

dx 


1 d_,dhs 
Pr  ^ Sc  Sc 


+ Q(x,t) 
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Combustion  Instability  Problem : Mode  Shapes 


First-Order  Upwind  Scheme 


TVD  Scheme;  Source  Treatment  TVD  Scheme;  No  Source  Treatment 
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^ 0^005 : — , — j 

0q r^2—  3 “ 5 6 7 8 

(d)  Controlled  variation  scheme;  Xr/H-4.43 

Streamfunction:  Backward-Step,  Turbulent  Flow 

t?o  - \V). '000 
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TVD  Schemes  can  be  Effectively  Extended  to  Sequential  Solvers 
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